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LIVIU I. NICOLAESCU 

ABSTRACT. We prove a Chern-Lashof type formula computing the expected number of critical points 
of smooth function on a smooth manifold M randomly chosen from a finite dimensional subspace 
V C C°° (M) equipped with a Gaussian probability measure. We then use this formula this formula to 
find the asymptotics of the expected number of critical points of a random linear combination of a large 
number eigenfunctions of the Laplacian on the round sphere, tori, or a products of two round spheres. 
In the case M = 5 1 we show that the number of critical points of a trigonometric polynomial of degree 
< v is a random variable Z v with expectation E(Z U ) ~ 2V0.6 v and variance var(Z„) ~ cv as 
v — > oo, c « 0.35. 
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Introduction 

Suppose that M is a compact, connected smooth manifold of dimension m. Given a finite dimen- 
sional vector space V C C°°(M) of dimension N we would like to know the average (expected) size 
of the critical set of a function v G V. For the applications we have in mind N ^> m. We will refer 
to V as the sample space and we will denote by fi(v) the number of critical points of the function 
v £ V. 
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More explicitly, we fix a Euclidean inner product h on V and we denote by S(V) the unit sphere 
in V. We define expected number of critical points of a random function in V to be the quantity 

M (M, V, h) := -L_ / /i(t;)|dS ft (t;)| 

(A.4) 1 /■ '-' 2 (Ato) 



/ e~^/x(v)|dl4(i;)|, V<r>0. 
Jv 



(27TC7 2 

In other words, //(M, V, /i) is the expectation of the random variable Zy,h 

S(V) Z Vth := M (t>), (Z) 

where 5(V) is equipped with the probability measure determined by the suitably rescaled area density 
determined by the metric h. 

Let us observe that we can cast the above setup in the framework of Gaussian random fields, [1,9]. 
Fix an orhonormal basis (^ a )i<a<N of (V, h), and a Gaussian probability measure on V, 



7 = — — ~e ^ \dv\ 
(2tt)t 



Then the functions £ a : V — > K, v h-> £ a ( v ) = («,^a), are independent, normally distributed 
random variables with mean and variance 1, and the equality 

V{x) = J2^(v)^ a (x) 

a 

defines an M-valued centered Gaussian random field on M with covariance kernel 

K v (x,y) = J^* a (aj)* a (y), Vx, y £ M. 

a 

Thus, we are seeking the expectation of the number of critical points of a sample function of this 
field. However, in this paper, most of the time, this point of view, will only stay in the background. A 
notable exception is Theorem 6.1 whose proof relies in an essential way on results from the theory of 
stationary gaussian processes. 

It is possible that all the functions in V have infinite critical sets, in which case the integrals in 
(/x ) are infinite. To avoid this problem we impose an ampleness condition on V. More precisely, we 
require that for any point x £ M, and any covector £ £ T*M there exists a function v £ V whose 
differential at x is £. As explained in [28, §1.2], this condition implies that almost all functions v £ V 
are Morse functions and thus have finite critical sets. 

The above ampleness condition can be given a different interpretation by introducing the evaluation 
map 

ev = ev v : M -»• V v := Hom(V, R), x i-> ev x , 
where for any x £ M the linear map ev^ : V — > R is given by 

evj.(«) = v(x), £ V. 

The ampleness condition is equivalent with the requirement that the evaluation map be an immersion. 

This places us in the setup considered by J. Milnor [24] and Chern-Lashof [9]. These authors 
investigated immersions of a compact manifold M in an Euclidean space E, and they computed the 
average number of critical points of the pullback to M of a random linear function on E. That is 
precisely our problem with E = V y since a function v £ V can be viewed canonically as a linear 
function on V v . 

The papers [9, 24] contain a philosophically satisfactory answer to our initial question. The ex- 
pected number of critical points is, up to a universal factor, the integral over M of a certain scalar 
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called the total curvature of the immersion and canonically determined by the second fundamental 
form of the immersion. 

Our interests are a bit more pedestrian since we are literally interested in estimating the expected 
number of critical points when dim V — > oo. In our applications, unlike the situation analyzed in 
[9, 24], the metric on M is not induced by a metric on V, but the other way around. We typically have 
a natural Riemann metric on M and then we use it to induce a metric on V, namely, the restriction 
of the L 2 -metric on C°°(M) defined by our Riemann metric on M. The first theoretical goal of this 
paper is to rewrite the results in [9, 24] in a computationally friendlier form. 

More precisely, we would like to describe a density \dp\ on M such that 

fj,(M, V, h) = i-— \ \dfj,\ = —f \dfi\, 

area(5(V)j J M <tn-iJm 

where \dp\ captures the infinitesimal behavior of the family of functions V . 

In Corollary 1.3 we describe such a density on M by relying on a standard trick in integral geom- 
etry. Our approach is different from the probabilistic method used in the proof of the closely related 
result, [8, Thm. 4.2], which is a higher dimensional version of a technique pioneered by M. Kac-and 
S.Rice, [1,21,32]. 

It is easier to explain Corollary 1.3 if we fix a metric g on M. The density \dp\ can be written as 
\dp\ = p g \dVg\, for some smooth function p g : M — > [0, oo). For x G M, the number p g {x) captures 
the average infinitesimal behavior of the family V at x. Here is the explicit description of p g (x) 

Denote by K x the subspace of V consisting of the functions that admit xasa critical point. Let 
S(K X ) denote the unit sphere in K x defined by the metric lion y. Any function v G K x has a 
well-defined Hessian at x, Kess x (v), that can be identified via the metric g with a symmetric linear 
operator 

Hess x (v,g) = Hess £C (f ,g) : T X M -»• T X M. 

We set 

A X (V):= f \ detRess x (v)\ \dS(v)\ [A = ] [ | detHes Sa! (t;)|e-W 2 |dv|. (A) 

The differential of the evaluation map at a; is a linear map A x : T X M — >■ V w , and we denote by 
J g (A x ) its Jacobian, i.e., the norm of the induced linear map K m A x : A m T x M -> A m V v . Then 



Pgi x ) 



J(At) 



and thus 



W 9 [X) 




_hf \ (Ml) 

|detHes Sg! H _ e J mKm \dV Kx (v)\ \dV g (x)\. 

K,. (27r) 2 / 

We want to emphasize that the density p g \ dV g \ is independent of the metric g, but it does depend on 
the metric h on V . In particular, the expectation p(M, V, h) does depend on the choice of metric h. 
For all the applications we have in mind, the metric h on V is obtained from a metric g on M in the 
fashion explained above. In this case we will use the notation p(M, V, g). Remark 2.2 contains a 
rather dramatic illustration of what happens when h is induced by a Sobolev metric other than I? . 
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The equality can be used in some instances to compute the variance of the number of critical 
points of a random function in V. More precisely, to a sample space V C C°°(M) we associate a 
sample space Va C C°°(M x M), the image of V via the diagonal injection 

A : C°°(M) -> C°°(M x M), Av(x, y) = v(x) + v(y), \/v G V, x,y G M. 

We equip Va with the metric Iia so that the map A : (V, h) — > (Va, ^a) is an isometry. Since 
/j,(A(v) ) = n(v) 2 we deduce that 

E ( Z V,h) = E { Z VA,h A )' 

where E denotes the expectation of a random variable, and are defined as in (Z). 

We also want to point out that if we remove the absolute value from the integrand Hess(u) in (fj-i), 
then we obtain a Gauss-Bonnet type theorem 

X (M) = 4ml T77TT ( I det Hess^) e"H 2 \dv\) \dV 9 (x)\, 

where x(M) denotes the Euler characteristic of M. 

Most of our applications involve sequences of subspaces V n C C°°(M) such that dim V n — > oo, 
and we investigate the asymptotic behavior of fJ,(M, V n ,g) as n — > oo, where g is a metric on M. 
One difficulty in applying (/l^) comes from the definition (A) which involves integrals over spheres 
of arbitrarily large dimensions. There is a simple way of dealing with this issue when V is 2-jet 
ample, that is, for any x G M, and any 2-jet j x at x, there exists v G V whose 2-jet at x is j x . 

Denote Sym^, the space of selfadjoint linear operators (T x M,g) — > (T x M,g). In this case, the 
lineal - map Hess : K x — > Sym^ is onto. The pushforward by Hessa, of the Gaussian probability 
measure f x on K x 

e 2 

1 X = , dim*, WkM\, 

(2tt) 2 

is a (centered) Gaussian probability measure 7^ on Sym^; see [20, §16]. In particular, 7^ is 
uniquely determined by it covariance matrix. This is a symmetric, positive define linear operator 

eJ:Symf ^ Symf . 

We can then rewrite (/2 X ) as 

n(M,V,h) = (27r)"t I \ det H \\drfl(H)\) \dV g (x)\. (7) 

JM J g (Al) \JSym*? J 

This is very similar to the integral formula employed by Douglas-Shiffman-Zelditch, [13, 14], in their 
investigation of critical sets of random holomorphic sections of (ample) holomorphic line bundles. 

In concrete situations a more ad-hoc method may be more suitable. Suppose that for every x we 
can find a subspace L x C K x of dimension £(x), such that for any v G K x , v _L L x we have 
Hess £C (t») = 0. Noting that dimi^a, = dim V — dim M = N — m and dimL^: = N — m — £(x) 
we obtain 



J e-\ v \ 2 \detRess x (v)\\dV(v)\ = \J ^ e^H 2 \dV{u)\ j x 
f e~\ w] ?\detlless x {w)\\dV(w) 
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N-m-i(x) f f i 12 

= tt 2 / |detHess 3 ,(iu)|e-l ,I 'l \dV{w) 



IM(M,V,h)=1T-T / 7T2-^±L\dV g ( X )\. (» 2 ) 
IM Jg{A ] x ) 



=:A(i 3! ) 

Using (a) we can now rewrite (/l^) as 

"1 A(-L a 

Our first application of formula (/x 2 ) is in the case M = S" l! ~ 1 and g is the round metric of 
radius 1 on the (d — 1) sphere. The eigenvalues of the Laplacian on S"^ -1 are 

X n (d) = n(n + d-2), n = 0,1,2,.... 

For any nonnegative integer n, and any positive real number u, we set 

V M := ker(A d - X n (d)) V v (d) := e y n>d . 

In Theorem 2. 1 and Corollary 5.2 we show that for any d > 2 there exists a universal constant > 
such that 

f x(S d - 1 , V v (d),g d ) ~ ^dimV^) ~ (d 2 ^ 1) , ^~ 1 as i/ -> oo. (A) 

The proportionality constant has an explicit description as an integral over the Gaussian ensemble 
of real, symmetric (d — 1) x (d — l)-matrices. 

In Theorem 2.3 we concentrate on the space V n) 2 of spherical harmonics of degree n on the 2- 
sphere and we show that as n — > oo 

If we denote by £ n the expected number £ n of nodal domains of a random spherical harmonic of 
degree n then, according to the recent work of Nazarov and Sodin, [26] , there exists a positive constant 
a such that 

Cn ~ ar ? as n — )■ oo. 

The estimate (B) implies that a < m 0.5773. The classical estimates of Pleijel, [31], and Peetre, 

[29], imply that a < 4 ~ 0.6916, where jo denotes the first positive zero of the Bessel function Jo- 

We next consider various spaces of trigonometric polynomials on an L-dimensional torus T L . To 
a finite subset M C Z L we associate the space V(M) of trigonometric polynomials on T L spanned 
by the "monomials" 

cos (mi #i H \-rriL0L), sin(mi0i H \-m L 9 L ), (mi, . . . , m L ) G M, 

and in Theorem 3.1 we give a formula for the expected number /i(M) of critical points of a trigono- 
metric polynomial in V(M). We consider the special case when 

M = := { (mi, . . . ,m L ) G Z L ; |m,| < i/, Vi = l,...,L} 

and in Theorem 3.2 we show that as v — > oo we have 

L 



KK)^(q) (IdetX^ xdimF(M^). 
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Above, (| detX|)oo denotes the expected value of the absolute value of random symmetric L x L 
matrix, where the space Sym L of such matrices is equipped with a certain gaussian probability 
measure that we describe explicitly. In particular, when L = 1, we have 

KK) ~ dim V(Mi) = 2vyjl, (E) 

while for L = 2 we have 

H(M 2 U ) ~ z 2 dimV(Ml). 

The proportionality constant z 2 can be given an explicit, albeit complicated description in terms of 
elliptic functions. In particular, 

z 2 0.4717.... 

In the case L = 1 we were able to prove a bit more. We denote by Z v the number of critical 
points of a random trigonometric polynomial in . Then Z v is a random variable with expectation 
E(Z U ) satisfying the asymptotic behavior (E). In Theorem 6.1 we prove that its variance satisfies 
the asymptotic behavior 

var(Z u ) ~ doou, (V) 

where 5^ is a positive constant « 0.35) described explicitly by an integral formula, (6.1). 

We also compute the average number of critical points of a real trigonometric polynomial in two 
variables of the form 

{acosx + bsinx + ccosy + dsiny + pcos(x + y) + qsin{x +y) }. 

This family of trigonometric polynomials was investigated by V.I. Arnold in [6] where he proves 
that a typical polynomial of this form has at most 8 critical points. In Theorem 3.6 we prove that 
the average number of critical points of a trigonometric polynomial in this family is 4^ ps 4.188. 
Note that the minimum number of critical points of Morse function on the 2-torus is 4, and the above 
average is very close to this minimal number. 

We then consider products of spheres 5' dl ~ 1 x S d2 ~ l equipped with the product of the round 
metrics x g d , 2 . In Theorems 4.4 and 5.6 we show that, for any di, d 2 > 2, there exists a constant 
-fQi,d 2 > sucn tnat ' f° r an y r — 1, as f — >• oo, we have 

xS d ^ 1 ,V u r(d 1 )®V u (d 2 )) ~K dli ^(dimV u r(d 1 )®V v ((k))"' {dl ' d2 ' r) , (C) 

where 

'1, (di-2)(da-2) = 0. 

(w) 



w(d 1 ,d 2 ,r) 



&Z%lt-l (di-2)(*-2)^0. 



Let us point out that for di, > 2, the function r i-)- ro(di, d 2 ,r), r > 1, is decreasing, nonnegative, 

lim tjj(di, d2, r) = = n(di,d 2 ) and ro(di, d2, r = 1) = 1. 

r->oo d 2 — 1 

In particular, 

w(di,d 2 ,r) < 1, Vr > 1. 

More surprisingly, 

w(di,d 2 ,r) = w(d 2 ,d 1 ,r) = 1 if (d\ - 2)(d 2 - 2) = 0, 

but this symmetry is lost if (di — 2)(d 2 — 2) 7^ 0. 

We find the asymmetry displayed in (C) + (117) very surprising and we would like to comment a 
bit on this aspect. 
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Observe that the union of the increasing family of subspaces W v>r = V u r(d\) (g> V^efe) is 
dense in the Frechet topology of C°°(M), M = S"* 1-1 x S^ 2 " 1 . The space C°°(M) carries a 
natural stratification, where the various strata encode various types of degeneracies of the critical sets 
of functions on M. The top strata are filled by (stable) Morse function. This stratification traces 
stratifications on each of the subspaces W v ,r and, as v — > oo, the combinatorics of the induced 
stratification on W ' vx captures more and more of the combinatorics of the stratification of C°°(M). 

The equality (C) shows that if r' > r > 1, the functions in W U:T - have, on average, relatively more 
critical points than the functions in W u y. This suggest that the subspace {W V)T ) captures more of 
the stratification of C°°(M) than W v y, and in this sense it is a more efficient approximation. The 
best approximation would be when r = 1, i.e., when the two factors S^ 1 participate in the process 
as equal spectral partners. Note that this asymmetric behavior is not present when one of the factors 
is S 1 . 

This heuristic discussion suggests the following concepts. Suppose that M is a compact, connected 
Riemann manifold of dimension m. Define an approximation regime on M to be a sequence of finite 
dimensional subspaces W, = (W v ) v >\ of C°°(M) such that 

Wi C W 2 C • • • 

and their union is dense in the Frechet topology of C°°(M). For any Riemann metric g on M, we 
define the upper/lower complexities of such a regime to be the quantities 

*/ w \ ,. logMM, W v ,g) . s. .. . log/i(M, W u ,g) 

k (W„g) :=hmsup— — - — — , n*(W;g) := liminf — — - — — . 

i/->oo log dim W u ;/->oo log dim W u 

Intuitively, the approximation regimes with high upper complexity offer better approximations of 
C°°(M). Finally, set 

k*(M) := sup K*(W.,g), k*(M):= inf K*(W.,g). 

W.,g w "9 

The above results imply that 

K*{S d ^), as*^ 1-1 x 5 d2_1 ) > 1, Vdi,d2>2, 

k^S* 1 ' 1 x S^- 1 ) < \fd!,d 2 > 3. 

a% — 1 

In particular, this shows that for any d > 3, we have 

k*(S 2 x S^" 1 ) = 0. 

In Example 5.4 we 1 construct an approximation regime (W n ) n >ion S 1 such that 

,. log y(S\W n ) 

hm — — — — = oo, 

n->oo log dim W n 

so that ^(S 1 ) = oo. 

Acknowledgements. I would like to thank Jesse Johnson for his careful proofreading of an earlier 
version the manuscript. 



The construction of the approximation regime in Example 5.4 was worked out during a very lively conversation with 
my colleague Richard Hind who was confident of its existence. 
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Notations 

(i) i : = ^=1. 

(ii) We will denote by <r n the "area" of the round n-dimensional sphere S n of radius 1, and by 
u n the "volume" of the unit ball in W 1 . These quantities are uniquely determined by the 
equalities (see [27, Ex. 9.1.11]) 

<r n _i = nu; n = 2— — y, r(l/2) = 0F, iff) 

where V is Euler's Gamma function. 

(iii) For any Euclidean space V, we denote by S(V) the unit sphere in V centered at the origin 
and by B(V) the unit ball in V centered at the origin. 

(iv) If Vo and V\ are two Euclidean spaces of dimensions no, n\ < oo and A : Vq — > V\ is 
a linear map, then the Jacobian of A is the nonnegative scalar J (A) defined as the norm of 
the linear map 

A k A : A k V -» A k Vi, k := min(no,ni). 
More concretely, if no < n\, and {e±, . . . , e no } is an orthonormal basis of Vq, then 

J(A) = (detG(A)) 1/2 , (J_) 
where G(A) is the no x ?iq Gramm matrix with entries 

Gij = ( Ae-i, Aej ) v ^. 

If ni > no then 

J {A) = J(A^) = (detC^t)) 172 , (J+) 

where A' denotes the adjoint (transpose) of A. Equivalently, if rfVoli G A n ' V* denotes the 
metric volume form on V\, and cTVoIa denotes the metric volume form on ker A, then J (A) 
is the positive number such that 

dVol = ±dVoU A A*dVoh . (J' + ) 

(v) For any nonnegative integer d, we denote by [a;]^ the degree d polynomial 

[x] d := x(x - 1) ■ ■ ■ (x - d + 1), 
and by Bd{x) the degree d Bernoulli polynomial defined by the generating series 

fptx fd 

d>0 

The d-th Bernoulli number is B d := B d (0), while the leading coefficient of B d (x) is equal 
to 1, and, 

B d+1 (u + l)-B d+1 



Tn d , Vz.GZ >0 . (5) 
d + l 

More generally, for any smooth function / : (0, oo) — > E and any positive integers u, m, we 
have the Euler-Maclaurin summation formula, (see [3, Thm D.2. 1] or [33, §7.21]), 

u-l fU m , 

n=l Jl k=l 

+ { ' ; / B m (x)f( m \x)dx, (EM) 
ml ./i 
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where bk denotes the /c-the Bernoulli number, bk ■= -Bfc(O), and B m denotes the associated 
periodic function 

B m (x) :=B m (x- [x\), VxGM. 
We will use one simple consequence of the Euler-Maclaurin summation formula. Suppose 
that f{x) is a rational function of the form 

/M = ^M, 

Pi(x) 

where Pq(x) and P\(x) are polynomials with leading coefficients 1 and of degrees do > d\. 
We further assume that / has no poles at nonnegative integers. Then 

t[ do - di + 1 

1. AN ABSTRACT RESULT 

Suppose that (M, 3) is compact, connected Riemann manifold of dimension m. We denote by 
I dV g I the induced volume density. 

Let V C C°°(M) be a vector subspace of finite dimension N. We set V y := Hom(V,R), and 
we fix a Euclidean metric /i = (— , — ) on V . We denote by S(V) the unit sphere in V with respect 
to this metric and by \dS\ the area density on S(V). The goal of this section is to give an integral 
geometric description of the quantity 

IM(M,V)=n(M,g,V 1 h):= [ m (v)\dS(v)\. 

area(S(V) ) J s{v) 

The significance of fi(M, V) is clear: it is the expected number of critical points of a random function 
v G S{V). 

To formulate our main result we need to introduce some notation. We form the trivial vector bundle 
V m '■= V x M. Observe that the dual bundle V M = V y x M is equipped with a canonical section 

ev:M->-V v , M9i4 ev x G V y , ev x (v) = v(x), G V. 

Using the metric identification V y — > V we can regard ev as a map M — > V. More explicitly, if 
{^a)i<a<N is an orthonormal basis of V, then 

a 

We have an adjunction morphism 

A : V X M -»• T*M, V X M 9 («, x) ^ := d^t; G T*M, 

where denotes the differential of the function v at the point x G M. We will assume that the 
vector space V satisfies the ampleness condition 

\/x G M the linear map V B v 1— ^ 4jjG T*M is surjective. (1.1) 

The assumption (1.1) is equivalent to the condition: 

the evaluation map ev : M — > V v is an immersion. (1.2) 

As explained in [28, §1.2], the condition (1.1) implies that for generic v G V, the restriction of the 
function v to K is a Morse function. We denote by fij\j(v) its number of critical points. 

For every x G M, we denote by K x the kernel of the map A x . The ampleness condition (1.1) 
implies that K x is a subspace of V of codimension m. Observe that the collection of spaces (K x ) x 
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is naturally organized as a codimension m-subbundle K — > M of V_ M , namely the kernel bundle of 
A. 

Consider the dual bundle morphism A' : TM — > V w x M. Using the metric identification 
V v — > V we can regard A^ as a bundle morphism A^ : TM — >• V_ M . Its range is K^ , the 
orthogonal complement to the kernel of A. Note that if {^ a }i<a<N is an orthonormal frame of V, 
x G M, and X G T X0 M, then 

N 

4 i = a(X4 fl )W)4 a ey. 

The trivial bundle V_ M is equipped with a trivial connection D. More precisely, we regard a section 
of u of V M as a smooth map u : M — > V. Then, for any vector field X on M, we define Dxu as 
the smooth function M — > V obtained by derivating u along X. Note that A * = D ev. 

We have an orthogonal direct sum decomposition V_ M = K L x K. For any section u of V_ M , we 
denote by u 1 - the component of u along K L , and by u° its component along liT . The shape operator 
of the subbundle K L is the bundle morphism S : TM ® .RT- 1 - — )• X defined by the equality 

S(X,w) := (D x u)°, VX G C°°(TM), w G C 00 ^). 

For every x G M, we denote by S x the induced linear map S x : T X M (g> K x — > K x . If we denote 
by Gr m (V) the Grassmannian of ?n-dimensional subspaces of V, then we have a Gauss map 

M3x^ S(x) := K x G Gr m (V). 

For a; G M, the shape operator S x can be viewed as a linear map 

S x : Ta-M -> Hom^,^) = T K ± Gr m (V), 

and, as such, it can be identified with the differential of 3 at x, [27, §9.1.2]. Any v G K x determines 
a bilinear map 

3 X ■ v : T X M (g> K x -)■ R, 3 X - u(e,u) = S x (e,u) ■ v, 

where, for simplicity, we have denoted by • the inner product in V. By choosing orthonormal bases 
(ej) in T X M and (uj ) of K x we can identify this bilinear form with an m x m-matrix. This matrix 
depends on the choices of bases, but the absolute value of its determinant is independent of these 
bases. It is thus an invariant of the pair (Ha,, v) that we will denote by | det 3^ ■ v\. 

Theorem 1.1. 

/i(M,V) = — !— / (/ \detS x -v\\dS(v)\)\dV g (x)\. (1.3) 

Proof. We denote by E x the intersection of K x with the sphere S( V) so that E x is a geodesic sphere 
in S(V) of dimension (N — m — 1). Now consider the incidence set 

E M := {(x,v) G M x S(V); A x v = } = {(a;, v) G M x S(V); v G E x }. 

We have natural (left/right) smooth projections 

M A- E M -A 5(F). 

The left projection A : .Eju — > M describes Em as the unit sphere bundle associated to the metric 
vector bundle Km- In particular, this shows that Em is a compact, smooth manifold of dimension 
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(N — 1). For generic v G S(V) the fiber p~ 1 (v) is finite and can be identified with the set of critical 
points of v : M — > R. We deduce 

F ) = rWu / ^ ^ \ dS ^ (L4) 

area (5(F)) / s(v) 

Denote by g E the metric on Em induced by the metric on M x S(V) and by \dV E \ the induced 
volume density. The area formula (see [16, §3.2] or [22, §5.1]) implies that 

f #p- 1 (w)\dS(v)\= f J p (x,v)\dV E (x,v)\, (1.5) 
Js(v) Je 

where the nonnegative function J p is the Jacobian of p defined by the equality 

P *\dS\=J p -\dV E \. 

To compute the integral in the right-hand side of (1.5) we need a more explicit description of the 
geometry of Em- 

Fix a local orthonormal frame (ei, . . . , e m ) of TM defined in a neighborhood IN" in M of a given 
point x G M. We denote by (e 1 , . . . , e m ) the dual co-frame of T*M. Set 

fi(x) := A\.ei{x) € V, i = 1, . . . , m, x G X. 

More explicitly, /j(a;) is defined by the equality 

(/i(x),») v = Se*«(as), Vi)£y. (1.6) 

Fix a neighborhood II C A~ 1 (Jsf) in M x S(V) of the point (a?o, Vo), and a local orthonormal frame 
u±(x, v), . . . ,un-i(x, v) over U of the bundle p*TS(V) ->• M X 5(V) such that the following 
hold. 

• The vectors iti(£C, u), . . . , u m (x, v) are independent of the variable i> and form an orthonor- 
mal basis of K^r. (E.g., we can obtain such vectors from the vectors fx(x), . . . , f m (x) via 
the Gramm-Schmidt process.) 

• For (x, v) G U, the space T V E X is spanned by the vectors u rn+ i(x, v), . . . , un-\(x, v). 

The collection tti (x) , . . . , u m (a;) is a collection of smooth sections of V_ M over X. For any a; G N 
and any e G T X M, we obtain the vectors (functions). 

D e m(x), D e u m (x) G V. 

Observe that 

E M nU= {(x,v) G It; l7 i (sB,«) = J Vi = l,...,m}, (1.7) 
where Ui is the function Ui : X x V — > R given by 

Ui(x,v) := (ui(x),v) v . 

Thus, the tangent space of E^f at (jc, w) consists of tangent vectors x®v G T^M © T V S(V) such 
that 

dUi(x, v) = 0, Vi = 1, . . . , m. 

We let a;;/ denote the m-form 

W£/ := dC/i A • • • A dU m G O m (U), 

and we denote by \\uju\\ its norm with respect to the product metric on M x S(V). Denote by \dV\ 
the volume density on M x S(V) induced by the product metric. The equality (1.7) implies that 

— 1 . 
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Hence 

j p \dV\ = 1 ^— [ \u u Ap*dS\. 
\\ u u\\ 

We deduce ^ 

Jp(x ,v ) = J p (x ,v )\dV\(ei, . . . ,e m ,ut, . . . ,un-i) 

= Ti—ti\ UJ u/\p*dS\(e 1 ,...,e m ,ui,...,u N -x) = -j— !— 7 \uj u (e l , . . . ,e m )\, r 
\\uuW \\uu\\ 'i&o,vol 

= :A u (x ,v ) 

Hence, 

/ #p-\w)\dS(v)\= [ -^\dV E (x,v)\. (1.8) 
JS(V) JE W^uW 

Lemma 1.2. We have the equality J\ = i, m , where J\ denotes the Jacobian of the projection 
A : E M -> M. 

Proof. Along It we have 

|dV"| = T1 rr- \oj v A cZVe I 

W^uW 

while (J^_) implies that 

\dV E \ = ^-\dV g AdS E J. 
Therefore, suffices to show that along It we have 

\dV\ = \ojuAdV g AdS E J, 

i.e., 

\u)u A dV g A dS Ex (e 1 , . . . , e m , tti, . . . ,ttjv-i) I = 1- 
Since dll^ujS) = 0, V/c > m + 1 we deduce that 

|w£/ A dT^ A dSs^iei, . . . , e m , u 1} . . . , wjv-l) I = • • • ,i*m)|- 

Thus, it suffices to show that 

\uu{u\, ■ ■ ■ ,u m )\ = 1. 
This follows from the elementary identities 

dUi(uj) = (ui,Uj) v = Sij, VI <i,j < m, 

where 5, L j is the Kronecker symbol. □ 

Using Lemma 1.2 in (1.8) and the co-area formula we deduce 

/ # P ~ 1 (w)\dS(v)\= [ ([ Au(x,v)\dS E M\) \dV g (x)\. (1.9) 
Js(V) Jm \Je x J 

=:J(a>) 

Observe that at a point (x,v) £ A _1 (X) C Em we have 

dUi(ej) = (D e .Ui(x),v) v . 
We can rewrite this in terms of the shape operator 3^ : T X M (g> K x — > K x . More precisely, 

dUi(ej) = (E x (ej,Ui),v) v . 

Hence, 

Au(x,v) = IdetHa, •« | , 
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We conclude that 

/ # P ~\v)\dS(v)\ = [ ([ \det3 x -v\\dS Ea (v)\) \dV M (x)\. 
Js(v) Jm KJEv J 

This proves (1.3) □ 

The story is not yet over. We want to rewrite the right-hand side of (1.3) in a more computationally 
friendly form, preferably in terms of differential-integral invariants of the evaluation map. The starting 
point is the observation that the left-hand side of (1.3) is plainly independent of the metric g on M. 
This raises the hope that if we judiciously choose the metric on M we can obtain a more manageable 
expression for [i(M, V). One choice presents itself. Namely, we choose the metric g on M uniquely 
determined by requiring that the bundle morphism 

A ] : (TM, g)^V x M 

is an isometric embedding. Equivalently, g is the pullback to M of the metric on V via the immersion 
ev : M — > V v = V. More concretely, for any x G M and any T X M, we have 

g x (X,Y) = (Alx,AiY) v . 

With this choice of metric, Theorem 1.1 is precisely the main theorem of Chern and Lashof, [9]. 

Fix x G M and a g-orthonormal frame(ei)i<j< m of TM defined in a neighborhood INT of x. Then 
the collection Uj = A^ej, 1 < j, is a local orthonormal frame of K L on N. The shape operator has 
the simple description 

a x( e i, u j) = ( D £i A^ ej ) . 
Fix an orthonormal basis (^ a )i<a<N of V so that every v G V has a decomposition 

v = ^2v a y a , v a G R. 

a 

Then, for any y G N, we have 

a a 

and 

( (D e .A^e j ) y ,v) v = ^v a {dl ie ^ a ) y = dl. e .v(y). 

a 

If v G K x , then the Hessian of v at x is a well-defined, symmetric bilinear form Hessa;(i>) : T X M x 
T X M — > R, i.e., an element of T*M ® T*M. Using the metric g we can identify it with a linear 
operator 

Ress x (v,g) :T X M^T X M. 

If we fix a g-orthonormal frame (ej) of T X M, then the operator Hessa;(v,5) is described by the 
symmetric m x m matrix with entries d\. e v{x). We deduce that 

|det 3^ • v\ = |det Hessa;(t7, g) | , G .Ba;. 
In particular, we deduce that 

M (M, V) = — !— /" (/ |detHess (B (t»^)||d5 a> (t;) > ) |d^ g (x)|. (1.10) 

Finally, we want to express (1.10) entirely in terms of the adjunction map A. For any x G M and any 
v G K x , we define the density 

Px,v ■ A m T x M -»• R, 
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p x , v {X x A • • • A X m ) = ^etid^x.vix))^ .^ | • (det( {rfX^Xfiv ) 1 < M < m )~ V2 

det(H«H ai («)(X i> X i )) la .< m | • (det(5(X 4 ,X J )) 1 ^ m )- 1/2 . 

Observe that for any g-orthonormal frame of T^M we have 

Po;,i;(ei A • • • A e m ) = | det Hess^, g) \ . 

If we integrate p x>v over t> G S(K X ), we obtain a density 

|d/i(£B, V)\ : A m T x M R 

that varies smoothly with x, and thus it defines a density \dp{— , V)| on M. We want to emphasize 
that this density depends on the metric on V but it is independent on any metric on M. We will refer 
to it as the density ofV. 

If we fix a different metric g on M, then we can express \dp(— , V)\ as a product 

|dM(aj,V)| = ^(aj)-|d^(aj)|, 

where p g = p g y : M — >• R is a smooth nonnegative function. 

To find a more useful description of p g , we choose local coordinates (x 1 , . . . , x m ) near x such that 
(d x i) is a g-orthonormal basis of T X M. Then 

P^(^ 1 A---A^ ro ) = |det(^. t ,( a; )) 1 ^ m | • (det((yita :Ei ,yita :Cj >) 1 < iJ .< rn ) _1/2 . 

Observe that the matrix ( d x . Xj v(x) ) 1<4 . <m describes the Hessian operator 

Hess x {v,g) :T X M^T X M 
induced by the Hessian of v at x and the metric g. 

The scalar ( det ( (A^d Xi , A^d Xj )v ) x<i j <m ) ^ i s precisely the Jacobian of the dual adjunction 
map A' x : T X M — >■ V defined in terms of the metric g on T X M and the metric on V. We denote it 
by J{Ax,g). We set 

Ax(V,g):= / | det Hess a .(u, g)\ \dS x (v)\. 

JS(K X ) 

Since 

\dV g (x)\(d xl A---Ad x J = l, 

we deduce 

p fl , v (aj) = & x (V,g) ■ J(A x ,g)-\ (1.11) 
We have thus proved the following result. 

Corollary 1.3. Suppose (M,g) is a compact, connected Riemann manifold and V C C°°(M) is a 
vector subspace of dimension N. Fix an Euclidean inner product h onV with norm \ — \h- Then 

(i(M,g,V,h) = tt^TYT / p(v)\dS h (v)\ = — [ ^ ,9 \ \dV g (x)\, (1.12) 

area (S(V) ) J\ v \ h=1 (tn-i Jm J(Al,g,h) 

where \dSh\ denotes the area density on the unit sphere {\v\h = 1}, and J(A x ,g,h) denotes the 
Jacobian of the dual adjunction map A x : T X M — >• V computed in terms of the metrics g on T X M 
and h on V. □ 



We will refer to the quantity p(M, g, V, h) as the expectation of the quadruple (M, g, V , h). 
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Remark 1.4. Let us observe that we have proved a little bit more. To every Morse function v G V 
we associate the measure 

fJ-v — ^ ^ &xi 
dv(x)=0 

where 6 X Denotes the Dirac measure concentrated at x. For every continuous function / : M — > M 
we set 



f_i(v,f):= fdfi v = V f(x), 
JM _, 77 n 



<fo(x)=0 

and we denote by /) the expection of the random variable S(V) B v \— > fi(v, /), 

Kv,f) = ^TT / Kv,f)\dS(v)\. 

area(S(V) J s(v) 

Arguing exactly as in the proof of Corollary 1.3 we deduce that, for any Riemann metric g on M we 
have 

= rWTT / v(v)\dS h {v)\ = — [ f*y ,9 \ f(x)\dV g {x)\. (1.13) 

axea{S{V)) J\ v \ h =i &N-\ Jm J(Ai,g,h) 

The resulting density on M 

-\dV g (x)\ 



a N - 1 J(Al s ,g, h) 
[ e-M 2 \detRess x {u)\\dV K „(u)) \dV g (z 



7r^2-J(Ax,g,h) 

is called the expected density of critical points of a function in V. As explained in the introduction, if 
V is 2-jet ample, then the above Gaussian integral over K x can be reduced to a Gaussian integral over 
Sym(T cc M). In this case, the resulting formula is a special case of [8, Thm.4.2] that was obtained 
by a different approach, more probabilistic in nature. □ 



Remark 1.5 (A Gauss-Bonnet type formula). With a little care, the above arguments lead to a Gauss- 
Bonnet type theorem. More precisely, if we assume that M is oriented, then, under appropriate 
orientation conventions, the Morse inequalities imply that the degree of the map p : Em — > S(V) is 
equal to the Euler characteristic of M. If instead of working with densities, we work with forms, then 
we conclude that 

X (M) = / —— dV g (x), 

&N-1 Jm J(A x ,g, h) 

where 



Xx(V,g):= / detB.ess x (v, g)dS x (v) 



>S(K X ) 

When M is a submanifold of the Euclidean space V, and we identify V with V v C C°°(V), then 
the above argument yields the Gauss-Bonnet theorem for submanifolds of a Euclidean space. □ 



We say that a quadruple (M, g, V, h) as in Corollary 1.3 is homogeneous with respect to a compact 
Lie group G if the following hold. 

• The group G acts transitively and isometrically on M. 

• For any function v G V, and any g G G, the pullback g*v is also a function in V. 

• The metric h is invariant with respect to the induced right action of G on V by pullback. 
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For homogeneous quadruples formula (1.12) simplifies considerably because in this case the func- 
tion p g y is constant. We deduce that in this case we have 

u(m, g, v, h) = Ax °Z ,g) r, ■ vol s( M )' (L14) 

<T N -iJ(Ai,. ,g,h) 

where xq is an arbitrary point in M. 

Let us observe that to any triple (M, g, V), V C C°°(M), we can associate in a canonical fashion 
a quadruple (M, g,V,h g ), where h g is the inner product on V induced by the L 2 (M, \dV g \) inner 
product on C°°(M). The expectation of such a triple is, by definition, the expectation of the associ- 
ated quadruple. We will denote it by u(M, g, V). We say that a triple (M, g, V) is homogeneous if 
the associated quadruple is so. 

2. Random polynomials on spheres 
As is well known, the spectrum of the Laplacian on the unit sphere (S , g^) C M. d is 

{X n (d) = n(n + d- 2); n > 0}. 

We denote by ^n,d the eigenspace corresponding to the eigenvalue \(d). As indicated in Appendix 
B, the space ^ n ,d has dimension 

A 2n + d- 2 fn + d- 2\ n d ~ 2 

M in, a) = ~ 2- — as n — > oo, 

n + d-2 \ d-2 J (d-2)! 

and can be explicitly described as the space of restrictions to 5 d_1 of harmonic homogeneous poly- 
nomials of degree n on U. d . For any positive integer v, we set 

V 

V v = V v {d) :=0Vn,d- 

n=0 

The space V v {d) can be identified with the space of restrictions to S^ -1 of polynomials of degree 
< v in d variables. Note that 

2v d ~ l 

N u := dimV v (d) ~ — — — as v -)• oo. (2.1) 
(d- 1)! 

The resulting triple (S ' ,gd, V u ) is homogeneous, and we denote by ^(S"^ -1 , y v ) its expectation. 
The goal of this section is to describe the asymptotics of u(S d ~ 1 , V v (d)) as v — > oo in the case 
d > 3. The simpler case d = 2 will be analyzed separately in Corollary 5.2. 

Theorem 2.1. For any d > 3 there exists a positive constant K = that depends only on d such 
that 

p{S d ^,V v {d)) ^ K d d\raV v (d) as v -»• oo. (2.2) 

In particular, 

logp(S d -\V„(d)) ~logdimV 1/ (d) as v ^ oo. 



Proof. For simplicity, we will write V„ instead of V^d). We will rely on some classical facts about 
spherical harmonics surveyed in Appendix B. For any integer d > 2, we denote by H n ,d the canonical 
orthonormal basis of constructed by the inductive process outlined in Appendix B and described 
in more detail below. 
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According to Corollary 1.3, it suffices to describe the density of V u at the North Pole p$ = 
(0, 0, ... , 0, 1) G S ,rf-1 . Denote by K u (po) the subspace of V„ consisting of functions for which po 
is a critical point. Note that 

dim K u (p ) = dimV v - dimS^ 1 = N„ - (d - 1). 

Near po we use x' = (xi, . . . , x d -i) as local coordinates so that 

Xd = s/ (1 — \x'\ 2 = 1 — ^\x'\ 2 + higher order terms. (2.3) 

Note that, at po, the tangent vectors d Xl , . . . , d Xd l form an orthonormal frame of T Pa S d ~ l . 

For any function / G C°°{S d ~ 1 ), we denote by Hess(/) the Hessian of / atpo, i.e., the (d — 1) x 
[d — 1) symmetric matrix with entries 

Hess(/)y = d%. x .f(p ), l<i,j <d-l. 

We set 

:={l,...,M(j,d-l)}, 

and we parametrize the basis %,<z-i as 

^>j,d-i = { ijf,^; /5 € -Bj.d }, 
where is a homogeneous harmonic polynomial of degree j in the variables x' = (xi, . . . , x d _i). 
For any integers j, n, < j < n, and any /3 G Bj,d, we define Z n j$ G C°°(S' d ~ 1 ) by 

Z nj ^(aj) := C7 n , i ,^](x d )y^( :E / ), G S**" 1 , (2.4) 

where j denotes the j-th order derivative of the Legendre polynomial P n)d denned by (B.2), while 
the universal constant C n j >d is described in (B.4). Then, for fixed n, the collection of functions 

{Z n ^eC co (S d - 1 ); 0<j<n, /3 G B j)d } 
is the orthonormal basis 'B n ,d- Any t> G V v admits a decomposition 

v n 

V = Yl Yl Y V n,j,fi Z n,j,0, V ntjjJ 3 G R, 
n=0 3=0 /3eSj, d 

so that 

v n 

Hess(u) = ^ ^ ^ -y nji/ g Hess ( Z nj - )/3 ) . 

n=0 3=0 /3eB, >d 

From the description (2.4) we deduce that 

Hess(Z nj , ) = 0, Vj>3. (2.5) 
Next, we observe that when j = we have M(0, d — 1) = 1 and 23o,d_i consists of the constant 

— 1/2 

function cr d _ 2 ■ We deduce 

"Bcd-l = {°"d-2 2 }' = i 1 }' 

Hess(Z nA i) = C nfitd cr~^ 2 2 Hess( P n ,dOd) )■ 

Using the equalities 

^n,d(i) = -Pm(I) + " 1) + higher order terms, 

we deduce 

Pn,d( x d) = PnA 1 ) + KA l ){ z d ~ l ) + hi g her order terms 
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P n (l) n,d ^ ■ \ x'\ 2 + higher order terms 

(B.3) 1 l, / d-3\. , l2 ,., 
= 1 — — [n + \)\n-\ — J \x | + higher order terms. 



This shows that 



Hess( P n , d (x d ) )=-l(n + l)(n+ U-i, (2-6) 

where t d -i denotes the identity (d — 1) x (d — l)-matrix. Hence 

Z n ,o,d G K u (p ), (2.7) 

Hess(Z ni o,i) =-\tr-^C nAd {n + l)(n + ^ (2.8) 

Similarly 

which implies that 

Hess(Z njli/ 3) = 0, Vn, p G B x = {1, . . . , d - 1}. (2.9) 

For any /3 G B 2yd , we denote by ffg the Hessian of Yp{x') at x' = G We deduce that 

^n,2,/3 € K v (p ), (2.10) 

Hess(Z ni2i/3 ) = C n , 2)d P2(l)^. (2.11) 
Using (2.5), (2.9) and (2.11) we conclude that 

V V 

Hess(t>) = X i> n ,2,/3 Hess(Z n;2 ,/3 ) + v n< o,i Hess( Z nA i ) 

n=2/3eB 2 , d n=0 
v 

= E (E^c„, 2 , d pS(i))^ 

1 f v d — 3 \ 
~ 2 <Td -2 2 ( E u «,o,iCn,o,d(w + 1) ( n H — ) I • 

\n=0 / 

The last equality can be rewritten in a more convenient form as follows. Define 

1 v d 3 
a = ao(u) := - -jV~ d -2 E ^-O.d^ + 1)(™H ^— )-^n,o,i G ^> ( 2 - 12 ) 

n=0 

and for /3 G -E>2,d> set 

a p = a p (u) := £ C n „ 2i ^ 2 ](l)Z nj2 ^ G V,. (2.13) 

n=2 

We deduce that 



Hess(t>) = (u,ao)l d -i + E^' a ^)- f ^- 



Note that the vectors clq, ap are mutually orthogonal, and they span a vector space L u of dimension 

* = £(rf) = M(2,d- 1) + 1 (B = } 
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Moreover, the conditions (2.7) and (2.10) imply that L u C K u (jpo). Define 



eo := r^rflo, := -^-a^, /3eB 2d , (2.14) 
l a o| | a/3 1 



d-3\* 



r = r (u) = |a | 2 = \<r d l 2 C n,o,d( n + 1 f( n + 



^ a— A / j ft,u,av 1 -/ y - 1 2 

n=0 

v 

Tp = rp(v) = \ap\ 2 = £ C^ id P®(i) 2 , /3 G B 2)d . 

n=2 

Note that the collection {eo, eg, /3 G -B 2i d} is an orthonormal basis of L„. For any i> G i^ i/ (po)> 
we denote by v its orthogonal projection onto L v , and we set 

v :=(v,e ), vp:=(v,ep), 

We deduce that for any i> G K v (po), we have 

Hess(i>) = Hess('u) = r^ 2 votd-i + ^^r^vpHp. (2.15) 

P 

For i> G K v (jpq) we set 

:= det Hess(u) I . 

Note that Q v {v) is positively homogeneous of degree d — 1. Using Lemma A.l in the special case 

n = &\mK u (po) = N u — (d — 1), m = dimLj, = £, uq = N v — (d— 1) — £, 
we deduce 

A X (V„)= [ Q u {v)\dS{v)\ =tr Nv - t - d f (1 - Ivl 2 )^'^' 1 Q u (v)\dV(v)\. 

JS(Ku(po)) JB{L V ) 

Using Lemma A.2 we deduce 



r( l+d-l \ r( N v -i-d+l 

B(i„) 2T(-f) -M-M 



(1 - |t>|0 a Q,(«) = V 2 - 2 / Qu(v) \dS(v)\ 



/ Q,(«) |dS(«)| = 1 2 V (2.16) 

•/soften)) r Hf 



Using the equality (a), we conclude 

h(K v ( P0 )) " "~ x " r(^f) 

Next, we compute the Jacobian of the adjunction map A' at po- We use the coordinates x' near po- 
For i = 1, . . . , d — 1 we have 

^^EE E] d Xi( Z n,j,p(Po))Zn,j,l3 = Yl 9 Xi( Z n,l,p(Po)) Z n>jt p. 

n=0j=0P£B jid n=ip & B hd 

Using (B.6), we deduce that Bi d = {1, . . . , d — 1} and for any /3 G Si d we have 

= O'd-^C'l&d-lZp, Z n,l,P = cr d--3C'ifi,d-iC n ^ td Pn :d {x d )xp. 

We deduce that 

V 

A\ d Xi =*-% 2 C lfi ^ l Y,C n ,l,dP'nA l ) Z nXi- (2-17) 

n=l 



20 LIVIU I. NICOLAESCU 

This shows that the vectors Ap d Xi , i = 1, . . . , d— 1, are mutually orthogonal and they have identical 
lengths 

V 



We deduce that the Jacobian of A Po is 



71=1 



Jv=r(v) A 2 L . (2.19) 



The equalities (1.14), (2.16) and (2.19) now imply that 

<Td-l _ Tj_ 



Using (a) we can simplify this to 

fi(S d -\ V v ) = K ll l v . (2.20) 

(vrr(^)^r(f) 

To obtain the asymptotics of V u ) as f — > oo we need to understand the asymptotics of the 

quantities 

r(v), r (v), rp(v), (3 G B 2) d-i- 
To achieve this, note first that (B.3) and (B.4) imply that 

r 2 P' m 2 - 1 (2n + rf-2)([n + d-3] d - 3 ) 2 (n + l) 2 / d-3^ 2 

^n.l.dW.dW 2 d- 2r (rf^l)2 [n + d-2] d _i 4 V 2 

= 1 ^2d-i(n) (r) 

2 d-i r (^i)2 B^^n) ' 

where A2d-\ (x) (respectively Bd-i(x)) is a monic polynomial of degree (2d— 1) (respectively d— 1). 
Using (Soo), we deduce that 



u d+1 , as v — > oo, 



2 

~ - 7T-777rfr^— - - - v d+ \ as */ -+ oo, (2.21) 



=1 n,M n ' Ml J ~ 2<*-i(d + i)r(%±) 2 

so that 

r 2 

2«*-ir(^) 2 < r d _ 3 (d + i) 

where 

2 (d-l)(rf-3)! 

2d - 3r( (£^ )2 • 

Invoking (B.3) and (B.4) again we deduce that 

r 2 m 2 - 1 (2n + d-2)([n + d-3] d _ 3 ) 2 A/n + 2 

nAdl J 2 d " 2 r(^i) 2 [n + d-l] d+1 V 4 V 2 

_ 1 A 2 d+z{n) 
~ 2 d +^Y{^Y B d+1 (n) ' 



d-3 
n H 



2 

2/ 
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where ^2^+3(2;) (respectively Bd-i(x)) is amonic polynomial of degree (2d+3) (respectively d+1). 
Using (Soo) we deduce that 

rp{u) = jyi 24 P^ d {lf ~ 2rf+ 3 (rf + 3 )r( ^_i )2 ^ +3 ' » ""►«>■ (2.22) 
Using (B.4), we deduce 

C laM + lf („ + i^f. (2 " + ;-j> [ ;^ 2 - 31 - (» + l) 2 (n + ^ 2 

= 2 j-3r ( ^i)2 - 4d+2( " ) ' <ro) 

where A d+2 (x) denotes a monic polynomial of degree d + 2. Invoking (5^) again we deduce that as 
i/->oowe have 

— V C 2 n o d (n + l) 2 (n + — V ~ — ^ v d+3 . (2.23) 



n=0 



Define 



r = lim f ( d+3 V (^), r« = lim v ( d+3 V«(z/), r= lim v ^ d+1 ^r(u). 



The precise values of these constants can be read off (2.21)-(2.23). Denote by the Euclidean 
space of dimension t = L) with Euclidean coordinates uq, up, j3 G B2.1i, and we set 

A u (u) = r Q {v)y 2 u Q l d ^ + Y,rp{v) 1/2 upHp, A^u) = f 1/2 u t d -i + (2.24) 

/? /? 

We can now rewrite (2.20) as follows 

li{S d -\V v ) = y l! / \ det A u (u)\\dS(u)\. 

(r{v)ir) — r(f ) ^5(ioo) 

The estimates (2.22) and (2.23) show that as v — > 00, we have 

(d+3)(d-l) 

I det A v (u)\ ~ v 2 I det (it) I, 
uniformly with respect to ii £ S(Loo). We deduce that as ^ — > 00, we have 

li{S d -\V v )~ \_\ J v d ~ x / IdetAooHUdSHI. (2.25) 

(vrf) — r(f) ^|«|=i 

This proves (2.2) where 



□ 



(vrf) — T(|)(d- 1)! V|«|=i V 2 



Remark 2.2. We want to analyze what happens to the above expectation if we change the L 2 -metric 
product on V v {d) to a new Euclidean metric so that the resulting quadruple (S d ~ 1 ,g,V u (d),h) 
continues to be homogeneous with respect to the action of SO(d). 



22 



LIVIU I. NICOLAESCU 



To perform such changes we use the fact that each of the spaces "% n $ is an irreducible representation 
of SO(d). Any sequence w = (w n )n>o of positive real numbers determines a Euclidean metric || — 
on V v (d) as follows. If 

V 

j) = ^D„gy,((i), v n g y n>d , 

n=0 

then we set 

- 1 

II II 2 — II II 2 

\\ v \\w '■- l^i 7^2 \\ v n\\ L 2^ S d-iy 
n=0 n 

In the sequel, we will choose the weights w of the form 

Wpv = {u + 1)P, pel. (2.26) 

The corresponding metric || — ||^ p is (equivalent to) the metric of the Sobolev Hilbert space H~ p 
consisting of distributions with "derivatives up to order — p in I?" . 

The quadruple (S^ 1 , g, V v {d), \\ — \\ w ) is homogeneous and we denote by n v (S d ~ l , V u ) its 
expectation. The collection 

{w n Z n ^£C°°(S d - 1 ); 0<j<n<u, p G d } 

is an orthonormal basis of V v with respect to the inner product h w associated to || — ||«,. Any v G V v 
admits a decomposition 

v n 

v = z2z2 /2 v n,j,pW n Z n ,j,l3, Vn,j,j3 G R. 
n=0 j=0 ^B 3 - d 

The arguments in the proof of Theorem 2. 1 show that for v G V u we have 

V 

P&B 2 , d n=2 



\ a d-2 ( z2 W n V n,0,l C n,0,d( n + ( n + ~ o~ ) ) ' ^-l" 
\n=0 / 



In particular, if Hess(i;) / 0, then the North Pole is a critical point of v, i.e., v G 2£"„(po)- Define 

1-1/2^ „ d — 3 

2 

n=0 

and for /3 G -E>2,d set 



a = a (u,w) := --a d ^ 2 2 ^2w n C nfi<d (n + l)(n^ — )z nfi ,i G V„, (2.27) 



fl/3 = a p {y,w) := ^n^d^ 1 )^ 2 '/ 3 G V u . (2.28) 

We deduce that 



Hess(v) = («,o )ld_i + ^(^a^i^. 

Define 

e := | — ,a , e^g := , -ap, /3 G £ 2 ,d, (2.29) 

l a o| I a/3 1 

^0 = r Q {v,w) := |a | 2 , = rp(v, w) := \ap\ 2 , j3 G B % d- 
For any i; G K u (po), we set 

:= I detHess(v)|, v = (v,e ), vp = (v,ep), /3 G B 2 ,d, 
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and we deduce 

Hess(» = Hess(u) = rl /2 v td-i + ^ r^vpHp, (2.30) 

P 

/ Q v (v)\dS(v)\ = , Q v C(v))\dS(v)\, (2.31) 

where denotes the orthogonal projection of t> onto the space L;, spanned by eo, e^. Similarly, we 
have 

V 

•A-l d Xi = o-^ 2 Ci )0) tj_i ^ io n C n) i )tf P^ )d (l)Z n) i )i . (2.32) 

n=l 

This shows again that the vectors yipo^, i = 1, . . . ,d — 1, are mutually orthogonal and they have 
identical length 

V 

n=l 

We deduce that the Jacobian of A Po is 

J y = r(z/)^. (2.34) 

If the exponent p in (2.26) is nonnegative, then using (r), (rg), (ro) and the Euler-Maclaurin summa- 
tion formula (5oo)> we deduce as before that as v — > oo we have 

r(z/, to) ~ Kt^+^P, rp{v lW ) ~ K 2 r d+l+2p , r (u, w) ~ A' 3 i/ ci+3+2 P, (2.35) 

where above and in the sequel we will use the symbols i^i , K 2 , . . . , to denote positive constants that 
depend only on d and p. This shows that 

H P (S d -\V^)~ K 4 dimV„. 

If the exponent p in (2.26) is <C 0, then a similar argument shows that 

Lip(S d -\ Vfj) ~ A' 5 , as i/ -> oo. □ 



We conclude this section with a computation suggested by the recent results of Nazarov-Sodin, 
[26]. 

Theorem 2.3. We denote by y n the eigenspace corresponding to the eigenvalue \ n = n(n + 1) of 
the Laplacian on S 2 , and we set /x(y n ) := n(S 2 , y n ). Then 2 

2 

~ ~~/= n2 as n ~ * 00 • (2.36) 
v3 

Proof. The computation is very similar to the computations in Theorem 2. 1 , but much simpler. We 
continue to use the notations in the proof of that theorem. In particular, K n denotes the space of 
harmonic polynomials in y„ that admit the North Pole po = (0, 0, 1) G R 3 as a critical point. 
An orfhonormal basis of y n is given by the 

Zn,j,p, < j < n, (3 £ Bj := B jj3 . 
2 Let us point out that -j= w 1.154. 
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Any y G V n admits a decomposition 



We conclude as in the proof of Theorem 2.1 that 

Hess(y) = ^2 Hess(Z ni2 „3) + y ,i Hess(Z n , ,i). 

We have 



(2.37) 



Hess(Z n ,o, 



1 



1/2 



C n ,o,3^(« + 1)12 



1 



2*r 2 ( 27r ) 1/2 



n(n + l)i/n+ -1 2 . 



In this case, the basis 1?2 consists of two elements, B 2 = {1,2}, and we have 

Z n ,2,l = Cn^^Pn), (x 3 )Yi (x 1 ,X 2 ), ^ n ,2,2 = CnjjPfflx^Yzfci,!?), 

where 

Yi(x 1 ,x 2 ) = Ci{x\ - ^2(^1,^2) = C 2 x\x 2 , 
and the constants Ci , C 2 are found from the identities 

1 = Cl I (x\ - x 2 2 ) 2 \ds\ = C 2 2 [ x\xl\ds\. 
Js 1 Js 1 

[ x\x 2 2 \ds\ = \ ( s\n 2 {26)\de\ = ^-, and I (x\ - x\) 2 \ds\ = I cos 2 {26)\de\ = vr, 
Js 1 4 Js 1 4 Js 1 Js 1 



Note that 



is 1 
so that 



We set 



Ci = tt~ 1/2 , and C 2 



7T 



1/2- 



ill := Hess(yi) 



7T 



1/2 



1 
-1 



We deduce that 



Hess(y) = ^,1^,2,3^(1)^1 + 2/2,2^,2,3^ (!)#2 - y Q ,i^m C nfiM n + 

2(7^ 



H 2 := Hess(y 2 ) 



>(2), 



7T 



1 



1/2 



1 

1 



7rV2 

We set 



^,2,3^(1)12/2, 





" 1 




" 1 " 




-1 


+ 2/2,2 


1 



1 



2/0,1 



o(n) := Cn, 2 , 3 P%l(l) = U(n+ly n + 2] 4 ) 



1/2 



b(n) :-- 



2(2tt) 1 /2 
1 



1 



ra(n+l)(n + -) 1 2 . 



1/2 



4^2 



n 1/2 

n(ra+l) (n+ - J . (2.38) 



Note that 



a[n) ~ — , on ~ — as n — > 00. 
w 64 v 7 32 



To ease the presentation, we set 



and we deduce that 



Hess(y) 



Ul'=V2,U u 2 :=V2,2, « = 2/0,1, 
2 



7T 



1/2 



au\ — bu au 2 
au 2 —aui — bu 
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Denote by L the space spanned by Z n 2,i, ^1,2,2, ^n,o,i> ^ := dim L = 3. It is contained in K n , and 
if y _L L, then Hess(i/) = 0. For 

V = UiZ n> 2,l + U 2 Zn )2 ,2 + uZ nt o,l € L, 

we have 

I det Hess(y)| = — \b 2 u 2 — a 2 u\ — a 2 u 2 \ ~ — — \2u 2 — u\ — u\\, n — > 00. 

7T 107T 

Arguing as in (2.17), we deduce that 

■4A = (9 Xl Z nXi (0))Z nAti , i = 1,2. 

Recall that 

^n,l,i = C C ni i i3 P' 3(3:3)^, 1 = Cq x 2 \ds\ = ttCq, 

Js 1 

so that 

1 /2 

I4a^l = Wn,l, 3 P; 3 (l)^ (B = } ^ 2 ( J^^ ) x in(n + 1) 

h \ V2 

n H — )n(n + 1) 



(4vr) 1 /2 ^ V '2 
Hence 

. n(n + l)(n+|) n 3 
J(Ko) = 4^ ~ ^ as « 00. 

Putting together all of the above and invoking (/i 2 )> we deduce that 



7? 2 



^ ~ area(S 2 ) x vr"^ x i / e"^ +u ' +u ^ \2u 2 - {u\ + u 2 2 )\ \dud Ul du 2 

4 JR3 



= -4/7 I e- {u2+u ' +u ^\2u 2 - {u\ + u\)\ \dud Ul du 2 \. 
Using cylindrical coordinates (u, r,9), U\ = r cos 9, u 2 = r sin 6, we deduce 

r p2n roo roo 

/ e -(« 2 +«?+«2) \ 2u 2 - {u 2 + u 2 2 )\ \duduidu 2 \ = / / / e- ( " 2+r2) |2u 2 -r 2 \rdrdud9 
Jr 3 Jo J -00 Jo 



00 roo 



= 2vr / / e~( u +r ^\2u 2 - r 2 \rdrdu = 4vr / / e"^ +r >|2u 2 - r 2 \rdrdu. 
J -00 Jo Jo Jo 

" v ' 

=:/ 

Hence 

MGU ~ 4vr- 1 / 2 /n 2 . 

To proceed further, we use polar coordinates u = tcosip, r = tsintp, 0<(p<^,t>0 and we 
deduce 



I 



roo I rn/2 \ 

/ / \2 cos 2 (p — sin 2 ip\ sin ipdip \ e~ t t 4 dt 
' 1 3 cos 2 (/? — 1 1 sin i/^fii^ 



00 

i 2 + 4 



(use the substitutions s = t 2 , x = cos <p) 



A • Qj3* 2 -l|dx) =^r(5/2) 



3\/3 2>/3" 
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Hence 



2 2 
v3 



□ 



Remark 2.4. Let us observe that for n very large, a typical spherical harmonic y G y n is a Morse 
function on S" 2 and is a regular value. The nodal set {y = 0} is disjoint union of smoothly embedded 
circles. According to the classical theorem of Courant [10, §VI.6], the complement of the nodal set 
has at most n 2 connected components called nodal domains. We denote by D y the collection of nodal 
domain, and we set 

%) := #V y < n\ 5 n :=7T-^ [ ''S(y)\dV \y)\ = l^S(y) \dS(y)\. 

In [26], it is shown that there exists a positive constant a > such that 

S n ~ an 2 as n — > oo. 

Additionally, for large n, with high probability, 6(y) is close to an 2 (see [26] for a precise statement). 

Denote by p(y) the number of local minima and maxima of y, and by s(y) the number of saddle 
points. Then 

Mi/) = p(v) + s (y), p(y) - = x(s 2 ) = 2. 

This proves that 

p(y) = 2 (Mv) + 2). 

For every nodal region D, we denote by p(y, D) the number of local minima and maxima 3 of y on 
D. Note that p(y, D) > for any I? and thus the number p{y) = J2oev p(Vi D) can be viewed as 
a weighted count of nodal domains. We set 

P0U:=7T-^ / e~\y\ 2 p{y)\dV{y)\. 

Theorem 2.3 implies that 

p(V n ) ~ as n OO. 

v3 

Since <5(y) < p(y), this shows that a < A=. □ 



Remark 2.5. We can use Remark 1.5 as a simple test for the accuracy of the computations in Theorem 
2.3. As explained in Remark 1.5, the Euler characteristic of S 2 is described by a integral very similar 
to the one describing n{^ n ). More precisely, we should have 

±2 = ±x(S 2 ) = area(5 2 ) x tt"^ x — L x ^ e Hyl 2 det Hess(y) |dV(y)| . 



The term J can be computed as follows. 

J= 1/2 , -tt / eH^ 2 detHe SS (y)|dy(y)| 

■K^l^nin + l)(n + ^) 

A simple application of the maximum principle shows that on each nodal domain, all the local extrema of y are of the 
same type: either all local minima or all local maxima. Thus p(y, D) can be visualized as the number of peaks of \y\ on 
D. 
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64 f r„,2. J _„2_ 1 _.,2i 



7r3/2 n ( n + i)( n+ 
(2.38) 1 



^ 2 + u t+<) ( 6( n )V - a{n) 2 {u\ + u 2 2 )) \dudmdu? 
{u*+ulW 2 )( 2 p(n)u 2 -a{n){ul+ul)) \dud Ul du 2 \, 



n 3/2 
where 

/3(n) = n(n + 1), a(n) = (n + 2)(n - 1). 
Arguing exactly as in the proof of Theorem 2.3, we deduce 

e _ (u 2 +u 2 +u 2 ) ^ _ a ( n )( n 2 + U2 j 2 | dudni(iu2 | 

noo 
e- (n2+r2) (2/3(n)n 2 - a(n)r 2 )r|drdu| 

(u = t cos ip,r = t sin </?) 

= 4ir / e _i / e _i t 2 (2/3(n) cos 2 — a(n) sin 2 sin <^|dtd<^| 



(x = cos (p, s = t ) 



oo />1 
.7 



ids x / ^ ( 2/3(ra) + a(n) )x 2 — a(n) ^jdx 



2 x 2 -(P(n)-a(n))=2^ 2 . 



2 3 

This confirms the prediction in Remark 1.5, namely, J = ±2. □ 

Remark 2.6. Most of the arguments in the proof of Theorem 2.3 work with minor changes for 
spherical harmonics of an arbitrary number of variables and lead to the conclusion 

^n,d) ~ Zdn^ 1 as n ^ oo, 

but the constant is a bit more mysterious. Here are the details. 
If 

n 

j=o /3e% d 

then 

Hess(y) = ^2 2/2,/? Hess(Z ni2 ,/3) + y ,i Hess(Z nA i). 
From (2.8), we deduce 

1 d 3 

Hess(Z„ i0 ,i) = --(Td!{ 2 C n , ,d(n + l)(^n-\ — ) 

1 _i /2 /(2n + d-2)[n + d-3] d _ 3 \ V2 (n + l)(n + ^) 
= -2^-2 x ^ ^ J >< ^ l^-i 

d+2 

1 n 2 

~ ~ 1/2 ^TT^-i as n ^°°- 

2 2 °d-2 1 W 

:a(d) 
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As in (2.4), we have 

Zn,2,p{x) ■= C nJ4 P^ d (x d )Y 2 ^(x'), x = (x',x d ). 
If we denote by Hp the Hessian of Yg at x' = 0, we deduce 

Hess(Z ni2|/3 ) = C nAd P^ d (l)Hp. 

Using (rp), we deduce 

d+2 

1 n 2 

Hess(Z n>2 ,/3) ~ d Hp as n -> oo. 

2^, i I— J 

=:&(<*) 

Arguing as in the proof of (2.17) we deduce 

Ap^dxi = 9 Xi Z n! i ) i(Po)Z ri! i ) i, 

where 

Z n ,l,i = C n xdP' n d{xd) ■ C d Xi, / Cjxf \dS(x')\ = 1. 

Js d - 2 



Hence 
We have 

Using (B.3) we obtain 
Using (B.8) we deduce 

Hence, 

so that 



i S d- 

•Aj, d Xi = C d C ny i td P' nd (l)Z ny i ti . 

1 rf-4 

C n ,i.d r — — n 2 as n->oo. 

' 2^r(^i) 

p 'nA l ) ~ as n -> oo. 
C 2 - 2 



7T 2 



KcAJ Fi— 7^ — wa as n -> oo, 

2 — r(^ i ) 

^(4o) ~ (^J r( d f i )( d-i ) - 

Denote by Z, the subspace of ^ n d spanned by the orthonormal collection of spherical harmonics 

{ Zn,o,i, Zn,2,p, P G B 2 ,d}- 

It has dimension 

dimZ = N d _x := dim Sym^S"^ 1 ) = 
Using (/x 2 ). Corollary 1.3 and the above computations we deduce 



n d-i 





-(d-i) , 


m 









where for 

V = ViZn,0,1 + ^2 y/3 Z n,2,/3 G 



CRITICAL SETS OF RANDOM SMOOTH FUNCTIONS ON COMPACT MANIFOLDS 



29 



we have 

Aoo{y) = yia(d)t d -x + ^ypb{d)Hp. 

P 

We interpret as an isometry from L to the space Sym rf _ 1 = Sym(T Po S , g) such that the 
collection 

a(d)ld-i, b{d)H/s 

is an orthonormal basis of Sym d _ 1 . We denote by g a ^ this 0(d — l)-invariant metric on Sym d _ 1 
and by | dV a j, | the associated volume density. We deduce 

/ e -M a |det4x>(v)| \dV{y)\ = [ e" |A| ^| det(A)| \dV a , b (A)\. 

JL v/Sym d _ 1 

On Sym fc _ 1 we have a canonical 0(d — 1) metric | — |* denned by 

\A\ 2 = trA 2 . 

We denote by \dV* \ the associated volume density. From (C.3), we deduce 

1 2 4T(^) 

\dV a , b \= ld \dVJ, ^= H{d _ lYHhR)Nd ^ R =~^- 

From (C.2), we deduce 

\A\l b = atrA 2 + f3(trA) 2 , 

where 

1 2 d+1 7r^ 



a 



b 2 R 2 r ^d±3) 

r d-2 



d-l\{d-l)a 2 b 2 R 2 J (d-l)l (d-1) T(^) 



2 d vr— / 1 2 \ 2 d 7r— /l 4 



(d-i) ^r(4fi)(d-i) r(^)y (d-i)r(^±i) V2 (d + i; 

2 d - 1 7rV(d-7) d-7 



2 

We deduce 



(d 2 -l)r(^) 8(d-l 



-a. 



■'Sym d _ 1 ■'Sym d _ 1 



As explained in Appendix C, the last integral can be further simplified to 

j e - atvA2 -^ trA ^\detA\\dV,(A)\ 

r 2 d—l 

= Z d / e"^-^^-) 2 TT \ Xi \ ■ \xi - Xj\ \dV(x), 

^ f=l l<i<J<d-l 

" v ' 

-.Id 

where tr(x) := x\ + • • • + x d ~i , and Zd is a positive constant that can be determined explicitly. The 
integral Id seems difficult to evaluate. The trick used in [18] does not work when d > 7, since in that 
case /3 > 0. The asymptotics of Id as d — > oo are very intriguing. □ 
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3. Random trigonometric polynomials with given Newton polyhedron 

Fix a positive integer L and denote by T L the L-dimensional torus T L := W L / \2ttZ) l equipped 
with the induced flat metric. Let 9 = (Q\, . . . , 0£) denote the angular coordinates induced from the 
canonical Euclidean coordinates on R L . For any m £ l/" we set 

( 2 1 / 2 



p(fh) 



(2^72 > 



A„l := p(rh) cos ^ m^- , B A := p(fh) sin ^ 

V=i / V =1 

The lattice 7L L is equipped with the lexicographic order -<, where we define fn -< n if the first non 
zero element in the sequence n\ — m±, . . . , — is positive. We define 6^ to be the positive cone 

C L : = {m £Z L ; -< m}. 

The collection 

{A ff }u{A^; meCiJu^mGeL} 

is an orthonormal basis of L 2 (T L ). A finite set M C Cl is called symmetric if for any permutation ip 
of {1, ... , L} we have 

(mi,...,m L ) e M^^fm^i),...^^)) e MU-M. 

For example, the set {(2, —1), (1, —1), (1, —2)} C C2 is symmetric. 
For any finite set M C Cl we define 



the scalars 



and the vectors 



V(M) := span{A^, B^; m, n<EM}, 
2 x - 

ajk = a jk (M) = L mjm k , (3.1) 

ajk = Oj fe (M) = y 2^ m i m k A rh G V(M). (3.2) 

( 27r ) 2 m S M 

If M is symmetric then the scalars are independent of j and we denote their common value by 
a(M). Similarly, the scalars ajk, j 7^ fc, are independent of j / fc, and we denote their common 
value by b(M). 

Theorem 3.1. Suppose Mc Gl is a symmetric finite subset of cardinality N > L. We set 

a:=a(M), b := /3(M), % = %(M). 

77?erc tlw following hold, 

(a) The sample space V(M) is ample if and only ifa^b. 

(b) Suppose that 

the vectors ajj are linearly independent. (#) 
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Denote by Sym L the Euclidean space of symmetric Lx L matrices with orthonormal basis (Hij)i<i<j<L> 
where Hj k is the symmetric L x L matrix with nonzero entries only in at locations (j, k) and (k,j), 
and those entries are 1. We denote by (—,—) the resulting inner product on Sym L . Then 

KM) = rju A ml e-^ c ~ lx ^\detX\\dX\, (3.3) 

(27r)2UJ(a-fe)— ( a + (L- 1)6) 1 Js y m L 

where for X = Yli<j x ij H ij 

\dX\ = JJ dxij 

i<j 

and C : Sym^ — > Sym^ is the symmetric linear operator described in the orthonormal basis (Hij) 
by the matrix 

Cij-ki = (aij,a k i) = , . L mimjm k m e . 

^ 71 ' meM 

Proof. We will compute /x(M) via the identity (1.14). Observe first that V(M) is invariant under the 
action of T L on itself, and the induced action on V(M) is by isometries. Let p = (0, . . . , 0) G T d , 
and denote by K p the subspace of V(M) consisting of trigonometric polynomials that admit p as a 
critical point. Set dj := de 3 , fj '■= dj\ p , j = 1, . . . ,d. We have 

A pfj = (p3^rn{Q)%n + djBm{6)Bfh \ \q = q = ^ m #)^' 

m€M mgM 

We have 

2 

G jk ■■= {-A-lfj^lfk) = t^Xl Yl m i mk = a ^( M ) 

^ ' meM 

Since M is symmetric we deduce that A p A P is described by the symmetric L x L matrix Gl{cl, b) 
whose diagonal entries are all equal to a, and all the off-diagonal entries are equal to b. We denote by 
Ai,(a,b) its determinant. We deduce 4 

A L (a, 6) = (a - b)^ 1 (a + (L - l)b) , (3.4) 

so that the Jacobian of A P is 

J(Al) = A d (a,b) 1 / 2 . (3.5) 

Observe that V(M) is ample if and only if the Jacobian of the adjunction map A P is nonzero, i.e., if 
and only if Ai(a, b) ^ 0. This proves part (a). 



If 



then 



meM 



i 

22 



djd k v(p) = ^2 arhmjm k . 



(2vr 



2 meM 



We deduce that 



i 

22 



Hessp(v) = ^ ^2 a rh m j m k H jk (3.6) 



(2vr 



2 j<fc \meM 



4 If Cl denotes the L x L matrix with all entries 1, then Gz,(a, 6) = (a — 6)lz, + 6Cl. The matrix C_l has rank 1 and 
a single nonzero eigenvalue equal to L. This implies (3.4). 
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Using the notations (3.1) and (3.2) we can rewrite the equality (3.6) as 

Ress p (v) = y2( v ^ij) H ij- 

Using Lemma A.3 and the equality dim V(M) = 2N we deduce that 



/ |detHess(t>)| \dS(v)\ = — |- [ e"H 2 | det Hess(«)| \dV(v) 
Js(K v ) r(iVj J Kp 

f e _b 5-|detHess(i;)| \dV(v)\ 



T(N) 

[ |detHess(i>)| & d L Kp \dV(v)\. 



2 

r+1, " 



T(N) 



=:7(M) 

We performed all this yoga to observe that J(M) is an integral with respect to a Gaussian density 
over K p . Denote by 7 = 7^ the linear map 

7 : K p -»■ Sym L , v H> Hess p (i>) = y~](v, dij)Hjy 

Since the vectors dy are assumed to be linearly independent, the map 7 is surjective. Clearly, 
det Hess(i>) is constant along the fibers of 7. As is well known (see e.g. [20, §16]) the pushfowrad 
of a Gaussian measure via a surjective linear map is also a Gaussian measure. Thus the density 



,,2 



\dlu\ ■=% I 6 Z Kv \dV(v)\ I , 



(27T)- 



2 

,|2 



is a Gaussian density on Sym L . Since the density — e A ~^ K \dV(v) \ is centered, i.e., its expectation 

is trivial, we deduce that its pushforward by T is also centered. The Gaussian density | d^y/i | is thus 
determined by its covariance operator 

C = Cm ■ Sym L ->• Sym L 

described in the orthonormal basis (Hij) by the matrix 

Cij-kt = (3ij,a k i) = mimjm k m £ . 
The symmetry condition on M imposes many relations between these numbers. We deduce 



J(M) = tjt^- / e~^ c 1X 'V\detX\ \dX\, s(L) := dimSym L . 

(2vr)V(detC) 1 /2 JSym L 

Using (1.14) we deduce 

T L.i dim K v 

vol(T L ) 2 _ 2- +1 7r — a 

^ = a < urn x — wm — J ( M > 

L dim K p £ 

= vol ( TL ) x 2 " ?7r ^ J(M) = ^ J(M) 
This proves (3.3). □ 
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We will put the above theorem to work in several special cases. Let us observe that the assumption 
(#) is automatically satisfied if M contains the points 

(1,0,. ..,0), (1,1,0,... ,0). 

Indeed, the symmetry of M implies that all the functions cos(#j) and cos(#j + 9j), 1 < i ^ j < L, 
belong to V(M) and the hessians of these functions span the whole space of L x L matrices. 
Suppose now that M = := n C, where v is a (large) positive integer, and A„ is the cube 

A^:={mGZ L ; |m y | < u, Vj = 1, . . . ,L }. 

Let us observe that 

K L V = M L V U (— M^) U {0}. (3.7) 
Among other things, this proves that M„ is symmetric. We want to investigate the behavior of /i(M^) 
as v — > oo. To formulate our next result we need to introduce additional notation. 
Let us observe that we have an orthogonal decomposition 

Sym L = V L (B7)i, (3.8) 

where D l consits of diagonal matrices 

T>l = span{i/jj; 1 < i < L} 

and 

T>j; = spanji/^; 1 < i < j < L }. 
For any real numbers o, b we denote by G^(a, b) the L x L-matrix with entries 



9ij 

Theorem 3.2. Let 



a, i = j 

b, i + j- 



:= {m £ Cl; \rrij\ < u, Vl<j<L}. 

Then, as v — > oo we have 

L_ 

K^v) ~ (f ) 2 < I det X\ >^ dim V(M£), (3.9) 

where ( I det X\ the expectation of I det XI with respect to the centered gaussian probability 
measure on Sym L with covariance matrix that has the block description 

C OQ = G L ( ?,1 ) film 



v2 • 



with respect to the decomposition (3.8). 
Proof. Let us first compute 

«M = ^) = 4 E -? ( = } T^r E ™? = ^ E ^l^ 1 ! 

v ; meM„ V 7 meA£ |mi|<i/ 



2 2^ + l L - 1 -A l9 2 2i/ + l L ~ 1 , , 1 T , 9 

y , ' > k 2 = y , ' B 3 (u + 1 ~ — r u L+2 as v -> oo. 

(2vr) L ^ 3(2vr) L dV 7 3vr L 

Similarly, we have 



b{y) = b(M%) = -^-t E rnim 2 . 
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The last sum is due to the invariance of with respect to the reflection 

(mi, 7712, . . . ,m L ) < — > (-mi, m 2 , . . . ,m L ). 

Thus, in this case 

A L (a,b) = a{u) L ~ — —^u L( - L+ ^ as v -> oo. (3.10) 
3^71"^ 

To compute the covariance operator C we observe first that, in view of the symmetry of M„ it suffices 
to compute only the entries 

Cn ;i j, i<j and C 12 -ij, i < j. 

We have 

1 4 2|A^| A 4 2(2i/ + 1) £ - 1 p , ... 1 L+ 4 

C/11-11 = i rr > 777 1 = — ; — > re = ; — -f BrAu + 1 ~ T v . 

11,11 (2ir) ^ 1 (2tt) l ^ 5(27r) jV ; 5vr L 

meAL k=l 

For 7 > 1 we have 

1 1 IA L -2| 

^ 1 \ " 2 2 1 \ " 2 2 il K \ " 2 2 

= 72^ ^ mi?ni = (2^ ^ mim2 = 72^ ^ mim2 

V 7 meAj' meA£ meAg 

Using the invariance of A v with respect to the reflections 

(mi, . . . ,rm, . . .,m L ) < — > (mi, ... , -m*, . . . ,m L ) (3.11) 
we deduce that for any i < j we have 

Cll,ij = 0. 

To summarize, we have shown that 



Xu = C u ,u = Cii ; ii ~ ^r-rv 1 ^ (3. 12a) 

= C«yj = Cnjj ~ ^T^ L+4 VI < i < j. (3.12b) 
C ii;j(t = 0, V», j</c. (3.12c) 

Next, we observe that 

^ 1 V- 2 2 4(2l/+l) L ~ 2 . 1 r, 4 

Ci2-i2 = 7 — rr > m m = v , ' B 3 (v + 1 = y u r u L+i . 

' (27r) ^ 1 2 9 2vr L tf 9vr L 

Using the reflections (3.1 1) we deduce that 

Ci2;y = 0, Vi<j, (i, j) ^ (1,2). 
With respect to the decomposition (3.8) the covariance operator has a bloc decomposition 



C 



9 3~ 
3~+ IK 



where 5F : —}T>l. The above computations show that 

g-=0, ! K = y v \^=y v \n^. 
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The operator S is described in the basis (Ha) of 2?l by the matrix Gl(x u , y v ) . We deduce that 

C = G L (x V: y v ) ®y v l(Ls = y u x (G L {z u ,l) © z„ = 



=:C V 

Using (3.12a) and (3.12b) we deduce that 

lim z v = -. 



detC~yP +L detG L f -,1 ) ~ f - ] ( - + L ) yP +L , as z/ -> oo. (3.13) 



We conclude that 

Using (3.3), (3.10) we deduce 

M(^) = — ^T/m,* / e-^< C " X ' X >|detX| 1^1 



(2vr)M2) l , H ^ :2i y|(2) + 2(detC ! ,) 1 /2 ^Sym 



1/2 

making the change in variables X = yj Y we deduce 



L L 



2 1 



(dim Symj +L) 



/i(M^) = 3 2 7r2 ^ 2 1/M 7 / e -l^" v > F )|dety| |cZF| 

^ v ul\ i(i+2) I( L )+^ , . /„ 1 II I 

(27T)2(.2) i y^2- a y2(2j+2( det C ^l/2 ^ym L 

= rnAP - / e-^^|dety| |dY|. 

(2vr)2(2)zy^-^(detC i ,) 1 /2 JSym L 

As f — > oo we have 

C, ^ Coo := Gl f -,1 ) 



Using (3.12b) we deduce that as v — > oo we have 

»(M V ) ~ Z Ll / L , Z L = — ^ / e -h(c'Jy,Y)\ detY \ \dY\. 

32(2vr)2(2)(detC 00 ) 1 /2 Jsym^ 

Since 

dimV(M^) ~ (2u) L as v — > oo, 

we deduce 

x dimSymL 1 — / e -l(^ )|dety||dF| _ 

V 6 7 (2tt) 5— 1 (det Coo ) 5 ^y"^ 

This proves (3.9). □ 

Let us apply the above result in the case L = 1. In this case consists of trigonometric polyno- 
mials of degree < v on S l , and Sym L = M. In this case we have 

/ii vi\ >/5 1 /" -sad. i, 1 f 00 _52£ , /6\* /I 

( det A W = — x —= I e is max = x —== / e « rcax = — \/ -. 



We deduce the following result. 
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Corollary 3.3. 



/i(M^,)~2W^, as u^oo. (3.14) 



When L = 2, the computations are a bit more complicated, but we can still be quite explicit. 
Corollary 3.4. 

^(M 2 ) ~ Z 2 dimV(My), Z 2 ~ 0.4717, as i/ -> oo. (3.15) 
Proof. We decompose the operators X G Sym 2 as 
so that det X = (xy — z 2 ). We write a := |, b := 1. Then 

11 1 

-(C^X,X) = -^^——(ax 2 + ay 2 - 2bxy) - -z\ det = (a 2 - b 2 ) 

1 f ~ 9.<J- h '^ (ax 2 +ay 2 -2bxy)-±z 2 2| 



" v ' 

=:I(a,b) 

As shown in Proposition A.4, the integral I(a, b) can be reduced to a 1-dimensional integral 

, / /-27T 2C 3 / 2 \ 

I(a,b) = V2n(a 2 -b 2 ) U - - d9 - 2va + 2vr j , 



where 
We deduce 

and 



c(0) := (a -b cos 20) 
2 - 2c 3 / 2 



(c + 2)V2 



d0-a + l« 1.7207... 



| x (|detX|) Uoo » 0.4717... 



□ 



3/2 

Remark 3.5. The antiderivative of ^ 2 )i/2 can De expressed in a rather complicated fashion in terms 
of elliptic integrals. □ 

Still in the case L = 2, suppose that 

M = {(1,0), (0,1), (1,1)}. (3.16) 
The space V(M) was investigated in great detail by V.I. Arnold, [4, 5, 6]. 
Theorem 3.6. If M is given by (3.16), then 

47T 

/i(M) = — ss 4.188. 
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Proof. We rely on Theorem 3.1, or rather its proof. In this case L = 2, dim V(M) = 6. The 
collection { A\ t o, Ai t i, Ao,i } is an orthonormal system, and we denote by L the vector space they 
span. Note that L C K p , and Hess(t>) = if v G L 1 - n K p . We have 



1 „ 2 1 1 „ 1 



nl6M 



Then 

2^r 2 

We decompose u G Las 



3 



a-&= — a + (L-l)b = ^, J(Al) = (a-b)—(a + (L-l)b) 



V3 
2tt 2 ' 



v = xA 10 + yA 01 + zAi i. 



and we have 



Hess p (i>) 



2 i/2 
"^t7 



x + z 2; 
-z y + z 



det Hess p ('u)| = — ^\xy + yz + zx\, \v\ 2 = x 2 + y 2 + z 2 . 



Using (/i 2 ) we deduce 



vol^xS 1 ) 1 
u(M) = ^ x — 



-(x 2 +y 2 +z 2 ) 



xy + yz + zx\ \dS(x,y,z))\ 



-(x 2 +y 2 +z 2 ) 



xy + yz + zx\\dxdydz\. 



The quadratic form Q(x,y,z) = xy + yz + zx can be diagonalized via an orthogonal change of 
coordinates. The matrix describing Q in the orthonormal coordinates x, y, z is the symmetric matrix 



1 1 

1 1 
1 1 



and its eigenvalues are 1, — |, — 5. Thus, for some Euclidean coordinates u,v,w,we have 



Q = -(2u 2 -?; 2 -u; 2 ) 



and therefore, 



The above integral can be computed using cylindrical coordinates (u, r, 0), 

v = rcos9, w = rsm9. 



r= ( v l +w 2)l/l 



We deduce 



27T /"OO /'OO 

/ / e -(u 2 +r 2 )\ 2u 2 - r 2 \rdrdudO = 2vr 

'0 J-ooJO 

(u = t cos tp, r = t sin y0<93<7r, t>0) 



00 />oo 



00 JO 



e" (n2+r2) |2-u 2 -r 2 |r(ir(ind 



/"OO / /*7T \ 

2tt / I / |2 cos 2 — sin 2 (p\ sin ^cfy? I e~ l t A dt 
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27T i [°° e -*y d A . 1 | 3x 2 _ 1]dx \ = n n°° e s s 3/2 d \ . 1 | 3x 2 _ 1]da : 



„. r(v2) .(vr).^. 



Hence, 

1 47T 

MM) = —/= y « 4.188. 



□ 



Remark 3.7. The typical trigonometric polynomial t G T(M) is a Morse function on S 1 x S 1 , and 
thus it has an even number of critical points. Morse inequalities imply that it must have at least 4 
critical points. We see that the expected number of critical points of a polynomial in V(M) is very 
close to this minimum, and that V(M) must contain Morse functions with at least 6 critical points. 
Arnold proved in [6] that the typical function in T(M) has at most 8 critical points. 

A later result of Arnold, [5, Thm. 1] states that a generic trigonometric polynomial in V(M) has 
at most 6 critical points. However, there is an elementary, but consequential error in the proof of this 
theorem. More precisely, a key concept in the proof is a (real) linear operator that associates to each 
holomorphic function / : C — > C a new function holomorphic function / defined by f(z) := f(z). 
Arnold states that if zq G C is a critical point of /, i.e., ^(zo) = 0, then it is also a critical point of 
/. Clearly this is true only if zq is real. For example, zq = i is a critical point of f(z) = (z — i) 2 , but 
it is not a critical point of f(z) = (z + i) 2 . □ 



4. A PRODUCT FORMULA 

Suppose that (M, g, V) is a homogeneous triple, m = dim M. We say that it is special if it admits 
a core, i.e., a quadruple (p,s, L, w), where p is a point in M, f = . . . , f m } is an orthonormal 
frame of T p M, L is a subspace of V and w £ V such that the following hold. 

Pi. The vectors A p f r , r = 1, . . . , m are mutually orthogonal. For any v G V we denote by 

Hess(i>) the Hessian of v atp computed using the frame /. 
P2. The subspace L is contained in K p = ker Ap, and for any v £ L 1 - we have Hess(i>) = 0. 
P 3 . w G K p n L L , w(p) / and \w\ = 1. We set L := L © span (w). 
P4. L 1 - C ker ev p . 

Remark 4.1. (a) The importance of a core stems from the fact that in applications we often have 

dim L <C dim V. 

We regard Hess as a linear map 

K p — > Sym(TpM) := symmetric linear maps T p M — > T p M. 

We observe that L D (ker Hess)- 1 -, so we would expect the dimension of L to be at least as big as 
( m ^~ 1 ) = dim Sym(TpM). In many applications, dim L is only slightly bigger than C^ 1 ) - 

(b) The conditions P3, P4 can be somewhat relaxed. We can define a core to be a subspace 
L C K p that contains ev p G V and satisfies P2. For example, if ev p G K p , we can choose L to be 
the sum between the line spanned by ev p and the orthogonal complement of ker Hess, but this space 
may be difficult to get a handle on in practice. For reasons having to do with the applications we have 
in mind, we prefer to work with the above more flexible definition. □ 
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Suppose that (p,s,L,w) is a core of the special triple (M,g, V). A basis of V is said to be 
adapted to the core if it can be represented as collection of functions Yj G V, j G J, where J is a set 
of cardinality dim V equipped with a partition 

J = {c}UlUl* UR m 

such that the following hold. 

• The collection (Yj)j & j is an orthonormal basis of V. 

• The collection (Yj)j G j is an orthonormal basis of L. 

• Y c = w. 

• The collection {Yj; j £ {c} U I U J* } is an orthonormal basis of K p . 

• R m = {1, . . . , m} and 

Y r = - — t — -yij/ r , Vr £ R m - 

\Apfr\ 

For such a basis, we write / := {c} U /, so that the collection (Yj)j e f is an orthonormal basis of 

L. 

Proposition 4.2. (a) Suppose that (M,g,V) is a special triple and (p, f , L,w) is a core of this 
triple. Set m : = dim M, i = dim L and N = dim V. Then 

vol (M)r( e+m ) f 

fjL(M,g,V)= l+ 9 : J \ 2 ' / \&A,H{v)\\dS{v)\ (4.1a) 

= ^ti M \ I e- H2 \^H{u)\\dV(u)\. (4.1b) 

(b) Suppose (M a , g a ,V a ), a = 1,2, are special triples with cores (p a , f a , L a ,w a ). Then the 
triple (Mi x M 2 , gi 052, Vi (8) V2) is special. The core is defined by the datum (p, /, L, w ), where 

p:=(pi,P2), / : =/ lu / 2 ) w^i,^) •"= xoi(xi)to 2 (x 2 ), 

erne? 

where L\ * L 2 denotes is the orthogonal complement ofw in L\ <g> L% Moreover 

I ev Pi,P2 I = I ev Pi I • I eV P2 I) ( 4 - 2 ) 

j <^> = • j«) • ( 'r e v ;''r. c ?:rir - <«> 

where J(S) denotes the Jacobian of a linear map between two Euclidean vector spaces. 



Proof, (a) Note that for any v £ K p , the Hessian H(v) of v at p depends only on the projection v of 
v on L. Using (A.2) we deduce 



N-m-e-2 . 



detH(v)\\dS(v)\=tr N - m -e-i / (1 - \x\ 2 )" "2 ' ~-\detH(x)\ \dV(x)\ 

S(K P ) JB(L) 

(A.3) cTN^^n^n^f^) r 

= ^tjh / |detff(aj)||d5(x)| 

2T( T ) Js(L) 

N-m-l „ . » ■ 

7T — 2 — r 

det dS^a?) . 



r(f) Js(L) 
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We deduce that 

KM,g,V) = ^#^-t- , / \det H(v)\\dS(v)\ 

= J^m.f |detfl"(t>)||AS(t>)|. 

2tt— nr=i i4/ri Js w 

(b) Choose an orthonormal basis (Yj*)j e j a of V Q adapted to the core (p a , f a , L a ,w a ), where 

J a = {c a }Ul a Ul* a UR ma . 

The collection 

{■^ii^Jo'i ,i2)eJix j 2 

is an orthonormal basis of <g> V2. 

Observe that HessfY"^ Y"?), the Hessian of Y}Yf at p = (pi,p 2 ), admits a block decomposition 

Yhp 2 )Re SS (Yl) A 



A* y/ i (p 1 )Hess(y j ? 



where for any (n, r2) € -R, the (n, r2) entry of the matrix ^4 is 

A ri r 2 = d f iY^(p x )d f 2 Y?(p 2 ), (n,r 2 ) G -R mi x i? m2 - 
The properties P2 and P4 show that if Hess (Y^Y?) 7^ 0, then 

either (ji,j 2 ) G A x I 2 \{(c 1 ,c 2 )} or (ji,j 2 ) G 4i x R 



7712 ■ 



This shows that if v _L L, then Hess(i?) = at (pi,p 2 ). The condition Z, C K Pl)P2 follows from the 
properties P 2 and P4 of special triples. 

Next observe that for any n G i? mi , we have 

KupiA = E E 

iieJi J26J2 
= E(%,^(^))^fe)Y r 1 1 ^. 

This proves that A^^fl-y {PI1P2) = (since Y^ (pi) = by P 4 ) and 

K^/n^l^^bOPEl^fe)! 2 



= l4 1 /n| 2 El^fe)| 2 = l4 1 AM ev P2 | 2 - 

J2=0 

We have an analogous formula for .Aj. P2 )/r 2 > r 2 e ^m 2 - The conclusions of part (b) of Proposition 
4.2 are now obvious. □ 



Example 4.3. An important example of special triple is (S V^), d > 3, where (5 1 

is the round sphere of radius 1, and V„ is the space spanned by the eigenf unctions of the Laplacian 

corresponding to eigenvalues X n = n(n + d — 2), n < v. A core can be constructed as follows. 
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As distinguished point, we choose the North Pole po = (0, . . . , 0, 1). Near po we use x' = 
(sci, . . . , Xd-i) as local coordinates, and we set 

f r = 9 Xr €= T po S 

— 1/2 

We choose w to be the constant function (r.'j . Finally, the subspace L is the space spanned by the 
functions e , denned by (2.12), (2.13) and (2.14), so that 

/ = {0} U B 24 . 

The properties Pi, P2 and P3 of a core are obvious. To prove property P4, we have to show that if a 
function v G V v is orthogonal to w, eo, ■ ■ ■ , ep, (3 £ -B2,d> then u(po) = 0. The function t; admits a 
decomposition 

v n 

V = EE E V n,j,P Z n,j,p- 
n=0 j=0 /3&B jtd 

Since Z n j : p(j?o) = if and only if j > 0, we deduce that 

V 

v (Po) = ^2v n>0 ,iZ nfi ,i(po)- 

n=0 

Note that w = Zo ; o,i- Since v _L w, we deduce 1*0,0,1 = 0, and therefore, 

n=l 

We now remark that (2.12) can be rewritten as 



1 v 

a = — r E Z n ,0,l(Po)Zn,0,l- 

71=1 

We deduce that 

V 

= -(t>,2a ) = ^*Vo,iZn,o,i(Po) = «(po)- 

n=l 

Note that (2.15) implies 

He SS (e ) = roH^ld-i, Hess( e/3 ) = r p (v) l l 2 H p . (4.4) 

From (2.17) we deduce 

V 

•Ay fr = & d-Z Clfi,d-1 E Cn,l,dP'n,d{^) Z nXr- ( 4 -5) 

n=l 

Using (2.18) and (2.21), we deduce that this triple has the additional property that 

\A\JA 2 = ■■■ = \Alj d J 2 = rid, v) ~ — -j — - — z/ +1 as 1/ -»■ 00. (4.6) 

N v ' 

=:r(d) 

Let us compute the length of the evaluation functional ev po = ev POjI , : V u — > BL We will use the 
notations in the proof of Theorem 2.1. We have 

v n v 

\^v,a 2 = Y.Y. E \ c n,j,dPnl^) Y j,m\ 2 = Y. E Ko^a^o) | 2 . 

»i=0i=0/9eB^ d n=0/3eB ,d 
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— 1/2 

In this case, the basis £o,d-i consists of a single constant function Yq = Y ^ = cr d _ 2 and 
p m-1 (B-6) (2n + d - 2)[n + d - 3] d _ 3 

Hence, 

v fn ,2 n 2 _ (2n + d-2)[n + d-3] d _ 3 

y o,^(uj c nAd (2-Ky- 1 ' 

and we conclude that 

1 - 2 
I ev ^ I" = J^yTT Y^ 2n + d-2)[n + d- 3] d _ 3 ~ J^yPi^ - 1) ^ as ^ °°' < 4 " 7 ) 



For r = 1, . . . , d — 1, we set 



The computations in the proof of Theorem 2.1 imply that 

d Xl U x (jp Q ) = d X2 U 2 ( Po ) = ... = d^Ud.^po) = c(d,u) = r(d,v) 1 ' 2 . (4.9a) 

d Xl U j ( Po ) = 0, = l,...,d-l, i + j. (4.9b) 

□ 



Theorem 4.4. Assume d\ , d 2 > 3 and fix r > 1. 77ze« f/zere exists a positive constant K that depends 
only on d\ and d 2 such that, 

zyVl, f2 — >■ oo and viv 2 r converges to a positive constant. The exponent w{d\,d 2 ,r) is described in 
(vj), 

2{dx - 3)r + 2d 2 + 2 



m(di,d 2 ,r) 



2{di - l)r + 2d 2 - 2' 



Proof. Choose cores (p a , f a , L a ,w a ) of (S 1 ^ 1 ,go,V l/a ) as indicated in Example 4.3. Next, 
choose bases adapted to these cores 

(Yj a ) jeJa3 J a = {c a }Ul a Ul* a UR da _ 1 , I a = {0} U B 2da , 

and set 

S= (Jj x J 2 ) \{(ci,c 2 )}, R = R dl ^xR d2 _ u I = SUR. 
Recall that for i a £ I a = {0} U B 2da , we have 



17' 



ep(v a ), i a = P£B 2}da 



where the functions eo, are defined by (2.14). Moreover, for r a € Rd a -i = {1, ■ ■ ■ ,d a — 1}, we 
have YJ? = C/ ra , where U r is defined by (4.8). 

We construct a core (p, f, L, w) of (S ' 1-1 x S^ 2-1 , V,^ <X> V V2 ) as in Proposition 4.2. Note that 
the collection 

{Y iuh :=Y % \Yl- (h,i 2 )el} 
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is an orthonormal basis of L. For v G V Ul ® V U2 we denote by Hess(-u) the Hessian matrix of v at 
p computed using the frame /. Note that if {i\,i2) G S, then 







Hess(T; lji2 ) 



^(p 2 )Hess(Y^) 

^(pOHess^) 
Using (4.9a) and (4.9b) we deduce that for (ri,r 2 ) G -R we have 

c(di,z/i)c(d 2 ,z/ 2 )A ri)r2 



r2; 



c(di,z/i)c(d 2 ,^ 2 )Aj ljr2 







c(di,vi)c(d 2l v 2 ) 



A 



Af j ,r2 



— :A ri , r2 

where A ri)r2 denotes the (di — 1) x (<i 2 — 1) matrix whose entry on the (ri, r 2 ) position is 1, while 
the other entries are 0. Thus, if 

V = ^ 1 ^11,12^1,12 ^ 
(11,12)6/ 

then 

Hess(^) = ^ v iu i 2 }less(Yi u i 2 ) + c(di,u 1 )c(d 2 ^2) ^ u riirj A riin . (4.10) 
(ii,»2)e5 (n,r 2 )eR 

To make further progress, we need to choose the basis (Y i a )i a ^i a of L a as indicated in Example 4.3. 
Using the notations in the proof of Theorem 2.1 we let 

I a = {0} U B 2Aa 

and the functions Yq respectively Yg, (3 G B 2 ^ a are equal to the functions eo and respectively 
defined by (2.14 ), (2.12), (2.13). More precisely, for a = 1, 2, we have 

1 ^ 



a 



o 



where 



and for (3 G -B 2 ,,i ct , we have 



where 



n=0 



dry — 3 



1 -1/2 



7 a r . 



a^ = Y,C n , 24 P^(l)Z n ^. 



n=2 



We set 



r (u a ) := |a-o | 2 , rp(y a ) := \a^\ 2 . 
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Using (2.22) and (2.23) we deduce that for any d > 2 there exist explicit positive constants ro(d), 
f i (d) such that 



ro(v a ) ~ r (d Q )^ Q+3 , as u a ->• oo 



^Oa) ~ n((i a )^ +3 , as ^ a 



oo. 



tig +3 



Let us observe that for any j3 G -E>2,d a we nave 

l^(p a ) = 0, Hess(F^) ( = 4 } rp {y a f^H p ~ n^vT^ Hp, as ^ Q -> oo. 
Moreover, 

d Q +3 

Hess(Y °) = r (^) 1/2 l dQ _i ~ ro(da) /2 f a 2 l^-i, as i/ Q ->■ oo. 

Next, 

Y a (p a ) = - 2 ,\ 1/2 <7dl-l £ C nAda n(n + ^-^)^n,o,i(Po), 



(4.11a) 
(4.11b) 

(4.12) 

(4.13) 



n=0 



where 



We deduce 



Y a (Pc 



7 , \ ( 2 - 4 ) n -1/2 

^n,0,l(PaJ - ^n,0,dc, a d a -2- 



2r (^) 1 /2 (TdQ _ 2 



2 rf - 1 r K)V2 <Tda _ 2 r(4 
1 



2 / n=0 



n=0 



cL — 3 



j- n( n + ) (2n + d a - 2)[n + d a - 3] 



rfa — 3 



Note that 



n=0 



2(4vr) — roK)V2^ 



dry — 3 



+ — (2n + ( i a - 2)[n + d a - 3] dQ _ 3 . 



^ n( n + ) (2n + d a - 2)[n + d a - 3] do 



d a + 1 



as v n -»■ oo. 



-i/rv as f„-> oo. 



Using (4.11a), we deduce that 

^ ~ "(47r)^ i (d a + l)r (d a )V2 " 
For i G we define the symmetric (d Q — 1) x (d Q — l)-matrix 

1 ' [Hp, i = peB 24a . 
Putting together all of the above, we deduce that for (ii,i 2 ) G S and v\, v 2 — > oo we have 

'0, i 2 G -B 2 ,d 2 or i! = ci 

z 2 = c 2 and i\ G 



(4.14) 



(4.15) 



r 2 (p 2 )Hess(^)~ < 



rfl+3 
Av \ 2 ^ 

^i^ 2 2 v 1 ' 2 H} i 2 = and n G Ii 



(4.16a) 
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n>i)Hess(^ 



fo, 

<%>+3 

CV 2 ff, 1 



«l G B 2 ,d! or « 2 = c 2 
i\ = c\ and i 2 £ J 2 



2 "22 

,A 2 fi 2 ^ 2 2 H i 2 > h = and i 2 G h 



(4.16b) 



where A,C are nonzero constants that depend only on d\ and d 2 , 1?^ is a nonzero constant that 
depends on i\ G Ji, , di, d 2 , and -Dj 2 is a nonzero constant that depends on i 2 G I 2 , d\, d 2 . 

Similarly, using the estimate (4.6) we deduce that for (ri,r 2 ) E i? and i/i, z/ 2 — >■ oo, there exists a 
nonzero constant i£ = E(d\, d 2 ) such that 



ril+l d 2 + l ^ 

Hess(Y rijr2 ) ~ Ev x 2 u 2 2 A ri)J . 2 . 



(4.17) 



Using (4.3), (4.7) and (4.6) we deduce that the Jacobian J(ui, v 2 ) of .At p v satisfies the asymptotic 
estimate 



J{v 1 ,v 2 ) = J(4). J(4J -|ev. 



|d 2 — 1 



(d 2 -l) 



-pi/ -\'-p 2 / I "'Pi I l cv Pll ^ "1 

where F is a positive constant. Assume now that 

Ui ~ t 2K1 , z/ 2 ~ t 2K2 , i -> oo, 
i.e., v\, u 2 go to infinity in such a fashion that 



(d 1 -l)(d 2 -l) 



1, r : = 



«i 



The assumption r > 1 implies that 
We have 



Kl > K 2 > 0. 



(4.18) 



d 2 -l rfj+3 

i/ 2 2 v x 2 

rf 2 -l di+3 

^1 2 v 2 2 

dl + l d 2 + l 
Z/i 2 IA, 2 



~ Ct Wn , wn = K\di + n 2 d 2 + 3ki - k 2 , 
~ C^ 22 , w 22 = Kidi + K 2 <i 2 + 3k 2 - ki, 



~ Ct W12 , CJi2 = Kidi + K 2 <i 2 + Kl + K 2 , 
1 2 ^2 



(d!-l)(d 2 -l) 
2 I/ 2 (>lZ/ 2 ) 2 ~ Ct 9 



where 



q = Ki(d\ - 1) + K 2 (^ - 1) + (ki + K 2 )(di - l){d 2 - 1) 
= Kid\ + K 2 d\ + (ki + K 2 )(did 2 — d\ — d 2 ). 

From (4.18) we deduce 

U 22 - Ull = — 4(«l - K 2 )<0, W12 - LUU = — 2(ki — K 2 )<0, 

P2 - = — (di + 3)ki + 4k 2 < -4(ki - k 2 ), 
pi - wn = -K 2 (d 2 - 1) < 0. 
Using (4.10), (4.15), (4.16a), (4.16b) and (4.17) we deduce that 

(B(V) + 0(1) ) < Wia_<,,U ^(Er 1 ,r a «n,r 2 A, 1)r2 + o(l) 



^a-am ^ ri , r2 A riir2 + o(l)) A^o^ 22 "" 11 ( l d2 _i + o(l) 



= :A(t) 
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where o(l) denotes a quantity that converges to as t — > oo, uniformly with respect to v G S(L), 
and 

B{v) = Y J B ll v llfl H} v 

h 

where are denned as in (4.16a). We set 

S(t) :- 

Observe that 

(B( V ) + 0(1) ) E(j2 ri r2 V rur2 A ritr2 + 0(i; 

^(Eri.ra u ri,ra^ri,ra + oi 1 )) A>«0,o( l^-i + o(l) ) 



f— 2(/ci— ks) 



= 5(t) 



We deduce that 



detff(w) ~t (dl+d2 - 2)wil |det5(t)| 2 -det 



r2 1Vi,r 2 ^ri,r 2 A)fO,old 2 -l 



= t (di-M a -2)wii-4d 2 (fci-fc 2 ) det j ^ oo ( u ). 

Let us point out that det H^v) is not identically zero. To see this, it suffices to choose v such that 
vq : q = 1 and all the other coordinates Vti,i 2 ^ trivial. In this case (4.15) and (4.16a) imply that 



Bold,-! 
A)l«fa-i 



It follows that 

^S^' 1 x 5 cZ2 - 1 ,F t 2 K1 ® V t 2 K2 ) ~ <#+^W-^ 2 (fci-fc 2 )-^ 
An elementary computation shows that 

(di + d 2 - 2)cjh - 4d 2 (fci - k 2 ) - q = 2d\Ki + 2d 2 K 2 - 6ki + 2k 2 . 
On the other hand, 

dim ® F t2K2 = dimV>«i x dimV t 2« 2 ( ~ 1} ^^2^+2^2-2^-2^ 
The desired conclusion follows by observing that 

2di«l + 2d 2 K 2 - 6«i + 2k 2 2(di - 3)r + 2d 2 + 2 . . « x 

; ; = —r~, ; ; = w(d-\ , a 2 , r), r= — . 

2«idi + 2K 2 d 2 - 2«i - 2k 2 2(di - l)r + 2<i 2 - 2 v ; k 2 



□ 
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5. Random polynomials on S 1 x S a \ d > 3 
For any m £ Z define X m : 5 1 — >• R by 

f (2^)~ 1 /2 ; m = o 

$ m (0) = } Tr-^cosfmfl), m<0 (5.1) 
I sin(m#), m > 0. 

The collection (^ m )rnez is an orthonormal Hilbert basis of L 2 (S l 7 d6). For any positive integer 
f > 0, we set 

T v :=span{$ m ; |m| <i/}. 
In other words, T v is the space of trigonometric polynomials of degree < v. 

Lemma 5.1. Let go denote the natural metric on S 1 of length 2tt. Then the triple (S 1 , go,T u ) is 
special. Moreover 

|ev P = !(, + !). (5.2) 

Proof. As base point we choose po = and the frame is / = {dg}. We denote by Kq the space 
of trigonometric polynomials that have as a critical point. Also, for any trigonometric polynomial 
t E T v we denote by Hess(i) the Hessian of t at 0, i.e., the 1 x 1 matrix Hess(i) := dgt(0)t. Note 
that 

A\d e = -tt~ 1/2 m $m. 

m>0 

In particular, 

/ v \V2 

J u := \A\de\ = vr" 1 / 2 HT m 2 ~ (Stt)" 1 / 2 ^ 3 / 2 as^oo. 



(5.3) 



\m.=l 



We set 

A simple computation shows that 
Next observe that 



\K d e\ 



9$p v (0) = J u - 



(5.4) 
(5.5) 



Hess(<3? r 



Thus, if 



then 



10, m > 

" j-m 2 $ m (0)l, m < 0. 

\m\<v 



m<0 



We now introduce 



so that 



Hess(t) = -k~ 1/2 (Y^ ^ 2 *m)l- 

i 

a = a u = -7T" 1 / 2 V] m 2 <I> r 

m<0 

Hess(t) = (t,a)l = |a„|(*, e„)l. 



e = e. 



a. 



(5.6) 
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We set 

1 " 1 1 

r v := \a u \ 2 = — N m 4 = — + 1) ~ — v b as i/ — > oo, 

7T ^— ' 57T 57T 

m=l 

and we observe that 

detHess(e„) = r 1 / 2 = (5vr)- 1/2 B 5 (v + 1) 1/2 ~ (5tt)- 1/2 i/ 5/2 as i/ -> oo. (5.7) 

We set c(v) = e u (0), and we observe that 

1 5 1 / 2 
c(„) = -— _ £ - 2 ^(°) ~ ~^T72 1/1/2 ^ ^ oo. (5.8) 

V ^ m<0 

We see that we can choose as core the quadruple (po,dg, L v , Yq), where L u denotes the 1 -dimensional 
space spanned by e„. The equality (5.2) follows from the identity ev = J2 m <o ®m(ty$m- □ 

Corollary 5.2. There exists a universal positive constant K such that 



H(S X ,T V ) ~ 2z^y|~ y|dimT^ as v 



oo. 



This agrees with our previous estimate (3.14). 
Proof. We use Proposition 4.2(a), and we have 



(use (5.3) and (5.7) ) 



/3 

t/-!/ as f — >• oo. 

□ 



Remark 5.3. Let us mention that, according to J. Dunnage, [15], the expected number of real zeros 
of a random trigonometric polynomial in 

T° = span{$ m ; < \m\ < v }. 

equipped with the L 2 -metric is ~ -%=v as v — > oo. The number of zeros of such a polynomial 
is a random variable ( v , and its asymptotic behavior as v — > oo has been recently investigated in 
great detail by A. Granville and I. Wigman, [19]. Let us observe the operator dg induces a linear 
isomorphism 

de : Tl -> T°. 

From this point of view we see that the expected number of critical points of a random trig polynomial 
in T° equipped with the L 2 -metric is equal to the expected number of critical zeros of a random 
trig polynomial in T°, equipped with the Sobolev norm \\u\\ H ~i. In Section 6 we will describe 
the asymptotic behavior as v — > oo of the variance of the number of critical points of a random 
trigonometric polynomial in T°. □ 

Example 5.4 (Approximation regimes with large upper complexity). Suppose 

^:{0,l,2,...,}-)-{0,l,2,...} 
is abijection such that (p(0) = 0. Define 

Tl := span{$ ± ^ (m) ; Q<m<v}. 
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We denote by /x(T^) the expected number of critical points of a random trigonometric polynomial in 
Tf}. A simple modification of the arguments used in the proofs of Lemma 5.1 shows that 

We want to construct 5 a permutation tp such that 

limsup — — — = oo. 

u^oo log dim Tl 

To do this we fix a very fast increasing sequence of positive integers (4)n>o such that 

4 = 0, = 2 » Vn > 1. 

mi 

For n > we set 

SVt := {4 + 1,4 + 2,..., 4+1 } . 

We consider the bijection cf> : {0,1,2,...,} — > {0,1, 2,...} uniquely determined by the following 
requirements. 

• 0(0) = 

• <f>(S n ) = S n , Vn > 0. 

• The restriction of ip to S n is strictly decreasing so that 

0(4 + 1) = 4+1, <j)(l n + 2) = 4+1 - 1 etc. 

We set 

v n :=£ n + l, W n :=Tt n . 
Note that the collection ( W n ) n >o is an approximation regime in the sense defined in the introduction, 
and 

dimW n = 24 + 3. 

We claim that 

lim i — T- — ?iTT =0 °- ( 5 - 10) 

n-s-oo log dim W n ) 

Indeed, for any positive integer k, we have 

^ 0(m) fc = ^ m fc + = P fc+1 (4) + i k n 2 \ 

m=l m=l 

where, according to (S), 

is a universal polynomial of degree k + 1. Using (5.9) we deduce 
Hence, 

lo gAt (W n )~2"log4, 

which proves the claim (5.10). □ 



'This construction was worked out during a lively conversation with my colleague Richard Hind. 
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Remark 5.5. Let us observe that for any positive v, the space 

T L V :=T V ® ■ ■ - T v C C°°{ S l x • • • x S\ ) 

s v y S v y 

L 

contains the space V(M£) of Section 3 as a codimension one subspace. The orhogonal complement 
of V(M£) in T£ is the 1 -dimensional space spanned by the constant functions. □ 

Theorem 5.6. For any d > 3 there exists a universal positive constant K = such that 
niS 1 x S^^pfdV^d)) ~ K d {dimT p ®V„(d)) as p,v oo. 



Proof. We will again rely on Proposition 4.2. We consider the core (0, dg, L p , w p ) of (Si,g s i , T p ) 
described in Lemma 5.1 and the core (po, f, L u , w u ) of (S 11 ' 1 , g S d-i , V „) described in Example 4.3. 
We form the core (p, f, L PiU , w) of (S 1 x S* ' -1 , 551 + g^d-i , T p <g) V^) following the prescriptions 
in the proof of Proposition 4.2. We have 

1 ~ p 
L p = span{e p }, e p = - — - a p = -vr" 1 / 2 ^ m 2 $ m (6>), ™ p (6>) = (2vr)- 1/2 . 

r^P 1 

r m=— 1 

As in the proof of Theorem 4.4, we choose a basis 

(Yj) j€ j, J a = {*}UlUl*UR d _ 1 , /={0}UB 2A , 
adapted to the core (po, f, L u , w u ). We have 



V 



' d-v 



For i G / = {0} U i?2,d Q . we have 

\e j g(i'), i = f3eB 2 , d , 
where the functions {eo(f)> e^(i/) are defined by (2.14). For r G -Rd-i = {1, . . . , d — 1} we have 



n=l 



We can now write down an orthonormal basis of L pv , 

{A iA :=w p Yi; iel}u{Z ■.= e p Y^}U{A i>2 :=e p Y i ; i e 1} U {B r := p p Y r ; r G 

where p p is given by (5.4). For any function v G V := T p W we denote by Hess(-u) its Hessian 
at (0,po)- We have 

" 



Hess^ii) 



(2^)- 1 / 2 Hess(F i ) 
Hess(A ii2 ) = 



Hess(Z) 



— 1/2 

<t Hess(e p ) 











Hess(e p )li(p ) 




Hess(p p y r ) 







p'J0)d x Y r (p )At 



e p (0)Hess(y 4 ) 
p' p (0)d Xr Y r ( Po )A 
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where A r denotes the 1 X (d — 1) matrix with a single nonzero entry equal to 1 in the r-th position, 
and Af denotes its transpose. In the sequel, the symbols Co, C\, . . . , will indicate positive constants 
that depend only on d. 

Using(5.5), (5.7) and (5.8) we deduce that as p — > oo, we have 

p' p (0) ~ CV /2 , Hess(e p ) ~ dp 5 , e„(0) ~ -C 2 p x l 2 . 
Using (2.22), (2.23) and (4.4) we deduce that as v — > oo, we have 



d+3 , , d+3 

Hess(Yb) ~ C % v~ \ d _ u Hess(y /3 ) ~ C 4 ^~H P (3 e B 2 , d . 



For i £ I we set 



Hi 



ld-i, i = 0, 

Hp, i = fie B 2 , d . 



We have 

Yp( Po ) = 0, V/3e5 2jd , 
while (4. 14) implies that as v — > oo, we have 

d-l 

*b(po) ~ -C 5 v 2 . 
Finally, using (2.18) and (4.9a), we deduce that as v — > oo we have 

d Xr Y r (po) ~ C 6 v^r. 
Putting together all of the above, we deduce that if 

V = Vi,lAi,l + ^ V i,2 A hl + zZ + ^ VrBr £ L P> U > 



iei 



iei 



then, as p, v — >■ oo, we have 
Hess(i>) ~ 



-C^ifi^V % + C 8 zp 5 CgEreflP 3 ^ '^rAr 

Cg E r6fi P 3/2 ^v r At v d ~¥ ( B ( W ) + pV 2 B x {v) ) 



where 



B (v) := C 10 v ,ild-l + ^2 C ii v i,i H i, Bi(v) := C" v , 2 ld-l ~ C 'u v h2 H i- 



Factoring out v ( 2 ' ) and then p 5 / 2 from the first row and the first column, we deduce that 

-C7V02 + C 8 z Co, YlreR P~ l w r /\ r 

C9 Erefl p- 1 ™^ u 2 (B (v) + p^B^v) ) 
(factor out v from the last (d — 1) rows and the last (d — 1) columns) 

-C 7 U 0) 2 + C 8 Z C 9 J2reR P~ ly r A r 

C 9 J2 reR p- 1 v r At B (v) + p 1 / 2 B 1 (v) 
(factor out p 1 / 4 from the last (d — 1) rows and the last (d — 1) columns) 

-C 7 V , 2 + C 8 Z Cg ~£reR P~ b/i VrA r 



d(d-l) 



det Hess(u) ~ u 2 p b det 



(d+4)(d-l) - d-l 

i/ 2 p " 1- 2 det 
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—C7V0 o + CgZ 

(d+4)(d-l) . I rf-1 , 

~ 1/ 2 p 0- *" 2 det 



To compute the Jacobian J p )iy of the adjunction map A po : 7(o )PO ) ( S" 1 x 5 d_1 ) — > T p <g> V u we use 
Proposition 4.2(a). We will denote by Kq, K±, . . . positive constants that depend only on d. 

According to (5.2) and (5.3), the Jacobian J p of A\ : TqS 1 — > T p , and the evaluation functional 
ev[J : T p — > R satisfy the p — > 00 asymptotics 

\ev p \ ~ K p 1/2 , J p ~ Ki/> 3/2 . 

Using (4.6) and (4.7) we deduce that the Jacobian J v of .Ap : T Po S d ~ l — ^ T p and the evaluation 
functional ev" = eVp Q : — >• R satisfy the ^ — > 00 asymptotics 

(d-i)(d+i) d_i 

,J U ~ A2f 2 , I ev I ~ K3V 2 . 

Using (4.3) ,we conclude that 

t t t I 0\d-l I 3 I (^H^ 2 ) 

</p,u = Jp • ^ • I ev | I eV I ~ ^4P 2 2^2 

We conclude that as p, v — > 00, we have 

□ 

6. The variance of the number of critical points of a random trigonometric 

polynomial 

The statistics of the zero set of a random trigonometric polynomial is equivalent with the statistics 
of the zero set of the gaussian field 



1 v 

(t) = y (a m cos mt + b m sin mt), 



'TTV 

m=l 



where a m ,b m are independent normally distributed random variables with mean and variance 1. 
This is a stationary gaussian process with covariance function 



1 " 

°M = — E 



cos mt. 

m=l 



The statistics of the critical set of a random sample function of the above process is identical to the 
statistics of the zero set of a random sample function of the stationary gaussian process 

j. , \ drjviit) 1 . , . 

qv\t) = — : = — > (— ma m sin mt + mo m cos mt). 

at v^v z — ' 

Equivalently, consider the gaussian process 

$„(t) = ^— _ (mcmcos + md m sin (J^jj 1 

where c m , d m are independent random variables with identical standard normal distribution, and the 
random variable 

Z v := the number of zeros of <& u (t) in the interval [—ttu, ttu\. 
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Note that the expectation of Z v is precisely the expected number of critical points of a random trigono- 
metric polynomial in the space 

span{cos(m#), sin(m#); m = 1, ...uj, 

equipped with the inner product 

(■2-K 

(u,v) = / u{6)v(9)d9. 
Jo 

The covariance function of is 

f mt s 

- y m~ cos 

7Ti/ c 



m=l 



\ V 



The Rice formula, [11, Eq. (10.3.1)], implies that the expectation of Z v is 



E(Z„) = 2u 



\ 2 {v)\ 2 



where 



V V 

AqH = Ru(0) = — j V m 2 and \ 2 {u) = -<(0) = —g V 



m 



This is in perfect agreement with our earlier computations. We let E((), and respectively var(C), 
denote the expectation, and respectively the variance, of a random variable (. The following is the 
main result of this section. 



Theorem 6.1. Set 



A := lim \ (u) = \, A := lim \ 2 {v) := -, 



Then for any t £ M the limit lim^oo R y {t) exists, it is equal to 

R 00 (t)--Lf r 2 C o S rdr = [ A 2 cos(\t)d\, \/t G R, 
t Jo Jo 



and 



where 



and 



i 2 r°° ( a \ 

lim - var{Z u ) = 5^:=- ( f^t) - J- ) dt + 2 ] 



(6.1) 



/«>(*) 



(^o ~~ ^0)^2 - Aq(-Rto) 2 



(Aq — i? 2 -,) 



- Plo+ Poo arcsin ) 



ro ° (A§ - i&)A 2 - A (i^) 2 ■ 
Moreover, the constant 6^ is positive. 6 

Proof. We follow a strategy inspired from [19]. The variance of Z v can be computed using the results 
in [1 1, §10.6]. We introduce the gaussian field 



*(ii,t 2 ) := 



^Numerical experiments indicate that <5oo ~ 0.35. 
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Its covariance matrix depends only on t = ti — t\. We have (compare with [19, Eq. (17)]) 

A R v {t) R' u (t) ' 

„ R u (t) A -R' v (t) T A B 

~ W " -R' u (t) A 2 -r'lit) "• [ St C 

_ R' u (t) -B»(t) A 2 

As explained in [30], to apply [11, §10.6] we only need that S(i) is nondegerate. This is established 
in the next result whose proof can be found in Appendix D. 



Lemma 6.2. The matrix S(i) is nonsingular if and only ift 2ttv'L. 
For any vector 



x 



%4 



G 



we set 



Pti,t 2 (x) :-- 



4 7 r 2 (detH) 1 /2 



g 2 ( a x i x ) 



Then, the results in [9, §10.6] show that 



As in [19] we have 



E{Zl)-E(Z u )=f ([ \ym\-Pt 1 ,t 2 (0,0,yi,y 2 )\dy 1 dy 2 

mr 1 : 



\dt1dt2 



* * 



,n = c -B^A~ l B. 



More explicitly, 











fi = 


c 




B^A^B 


A 2 




-Kit) - 












-Kit) 


A 2 












1 















A -J 


-R»it) 2 


#,(t) 











-ik(t) 




A 


2 — -R" 


iK? 




Ao Rv 




—R" A 2 






-Rl 


Ru Aq 



Kit) 
-Kit) 



1 - P 
-p 1 



where 



/ \2 



P- = P-v 



(Aq - R^)A 2 - A p^) 

Aq ~~ ^ 



fh 



Ki\ - R l) + iK) 2 Rv 



□ 



(6.2) 



(A2-i?2 )A2 _ Ao(jR /)2 • 

We want to emphasize, that in the above equalities the constants Ao and A 2 do depend on v, although 
we have not indicated this in our notation. 

Remark 6.3. The nondegeneracy of S implies that p u (t) 7^ and < 1, for all t 2-kvIa. □ 



We obtain as in [19, Eq. (24)] 

det S = det A ■ det U = p?{\l - i#)(l - p 2 ), ft -1 



1 



M(l - P 2 



1 P 
P 1 
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We can now rewrite the equality (6.2) as (t = t 2 — t\) 

y*+2„M + vl \ d y ld y 2 \\ \dhdt 2 \ 



E(Zl)-E{Z u )= / / \ ym \e ^(i-^) 



I v Kl v \J* 47T2 J ^ {X 2-R2 ){1 _ p 2 

yidy 2 \ \ fi\dtidt 2 \ 
4P~J V(Ag-i^)(l-^" 



/" f /" l I - xl+ ^ xl X 4 \dyidy 2 \ \ (Adtidtz 

/ / \ym\e w-r) 1 

Jl v xl„ \JR 2 



From [7, Eq. (A. 1)] we deduce that 

/ 2/12/2 e *(i-p*> 1 ^ ,f 21 = f- 1 + J__ arcsinp . 

Hence 

£(Z 2 ) - S(Z„) = / r i^Jl- p 2 + p&icsmp) \dtidt 2 \ 

Jl„xl„ (\%-R?,)2 V / 



r„ (A 2 - 

= 1 /■ (Ag - ig)A 2 - yfr) 2 (^p2- + paics - mp ) \ dtldt2l (6 .3) 

V v ' 

=:/*(*) 

The function f v {t) = f v (t 2 —t\)is doubly periodic with periods 2irv, 2irv and we conclude that 

2u f nu 

E([Z U } 2 ) := E{Z 2 V ) - E(Z U ) = — / f u (t)dt. (6.4) 

We conclude that 

«flr(Z„) = E{[Z V ] 2 ) + E{Z„) - {E(Z v )f = E([Z V ] 2 ) - [E{Z V )} 2 

2v r v ( „ , , A 9 \ , /A7 (6-5) 

To complete the proof of Theorem 6. 1 we need to investigate the integrand in (6.4). This requires a 
detailed understanding of the behavior of R u as v — > oo. It is useful to consider more general sums 
of the form 

a / \ 1 v "\ r TTlt . . 1 y > r . Vflt 

A ur (t) = — — > m cos — , B ur (t) = — — - > m sm — ,r>l. 
' K ' u r+1 ^ v ' u r+1 ^ v 

m=l m=l 

Note that if we set z := cos ^ + i sin We have 

1 JL 

m r z m . 



Observe that 

pWfri - 

We set 




(6.6) 



Aoo t r(t) = y r r cosrdr, B^^t) := ^ r r sinr(ir. 
Observe that it^ = ^2 and R^ = Ax>,2- We have the following result. 
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\A„, r (t) - A»,r(t)| + \B v>r (t) - Soo,r(*)| = 



max(l, t) 



, V< > 0, 



(6.7) 



where, above and in the sequel, the constant implied by the O-symbol is independent of t and v. In 
particular 



lim -L-C (i) = C r (t) : 



u^foo v 



r+1 



fr+1 



T r e tT dr, Vt > 0. 



(6.8) 



Proof. We have 



W = ^rE"'~(v)=FTt{(^)'~(^)} © 

m=l x 7 m=l v \ / 



:S„(t) 



The term S^i) is a Riemann sum corresponding to the integral 



and the subdivision 



f(r)dr, f(r) := r r cosT, 



<-<••• < ^ < t. 

v v 



of the interval [0, t\. A simple application of the mean value theorem implies that there exist points 

(m — l)t mt 



v v 



such that 



We deduce that 



rt " 

\ f(r)dr=J2 

J m =l 



f(0n 



111 



[ t f(T)dT-S„(t)= t -j2(Wm)-f(- 
V m=l ^ ^ 



Now set 



Observe that 



We deduce 



M(t) := sup \f'(T)\. 

0<T<t 



M(t) 



S u (t) - / f{r)dr 



0(t r ^), 0<t<l 
0(f), r > 1. 



< M(t) ■ -. 



This, proves the A-part of (6.7). The i?-part is completely similar. 



□ 



We need to refine the estimates (6.7). Recall that [m] r := m(m — 1) • • • (m — r + 1), r > 1. We 
will express C v%r (i) in terms of the sums 



d 



d r ( z — z 



m=l 



dz r If-' / df 1 



d r /l-z 



<fe r V 1-2 
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Using the classical formula 

^2S(r,k)[m] k , 



r 

r * 

m = 

k=l 



where S(r, k) are the Stirling numbers of the second kind, we deduce, 

c, r (o = E s ^ k ) D »A0 = E s ^ k ^ r ^k ( -737- ■ (6 - 9) 

k=l k=l 



Lemma 6.5. Set 9 := ±, and f(6) = Then 



(6.10) 



Proof. We have 



D ur (t) = z r Y ( r ) £(1 - - z)" 1 

E(-)" • 



_ /A „i/+l+r— j 



Using the identity 



we deduce 



. 2 . fa\ 
1 — e = - sm — e 2 
% V2> 



2sin(^) 

and 



(2sin0) r+1 6 6 f-T* Vi/ V2sin0yl 



Multiplying both sides of the above equality by we get (6.10). □ 

Lemma 6.5 coupled with the fact that the function f{9) is bounded on [0, |] yield the following 
estimate. 

t r+1 

—,D v Jt) = 0(l), Vi/, 0<t<nu. (6.11) 
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Using (6.11 )and the identity S(r, 1) = 1 in (6.9) we deduce that there exists K = K r > such that 
for any v > and any t E [0, 7rv] we have 



t r+l 



,r+l 



r-1 f+1 



3=0 



3=0 V 



so that 



-l r (c 1 ,, r (t)-z> I , >r (t)) 



< K r 



Using Lemma 6.5 we deduce 

t r+l 



lim 



v— >oo 1/ 



r+l 



ReA,, r (t) = I r {t) := i r r\ 2sin - • - e u ^ 



uniformly for i on compacts. The estimate (6. 12) implies that 



r+l 



lim ^— ReA,, r (t) =I r (t). 
We have the following crucial estimate whose proof can be found in Appendix D. 
Lemma 6.6. For every r > there exists C r > swc/j that for any u > we have 



1 



,r+l 



A/,r(t) " ^Ir(t) 



< 



C r t r+l - 1 

7r(t - 1) 



, V0 < f < TTV. 



V 



(6.12) 



(6.13) 



□ 



Using Lemma 6.6 in (6.12) we deduce 



v 



r+l "■>' v ' fr+1 ' v ' 



< 



C r t r+1 - 1 

~t r (t - 1)' 

r-t 



V0 < t < itv. 



(6.14a) 
(6.14b) 



I r (t) = t r+l C r (t) = / T r e tT dr. 

Jo 

Using (6.7) and (6. 14a) we deduce that for any nonnegative integer r there exists a positive constant 
K = K r > such that 

\C v Jt) - C r (t)\ < ^ f^"! - (6-15) 



Coupling the above estimates with (6.7) we deduce 



v t r (t- 1)' 



C vr (t) = C r (t) + O - , V0 < t < v 



(6.16) 



v. 



(6.17) 



where the constant implied by the symbol O depends on k, but it is independent of v. The last equality 
coupled with (6.6) implies that 

R?\t) = R&\t) + o(±Y V0<t< 

We deduce that, for any t > we have 

lim f v (t) = /«,(*) 

V— >oo 

where /j, is the function defined in (6.3), while 

(^o ~~ R"L)^2 - Ao(-R^) 2 



/oo(t) 



(*8 - i£ 



\A - P^o + Poo arcsin p c 



(6.18) 
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where 



A = A (oo) = lim AqH 



-, A 2 = lim \ 2 {v) = -, 



Poo(t) = lim p v (t) 



R'oci^O — Rqo) + (R'oo) R<x 

^oo- • (A2- J R2 o )A 2 -A ( J R^) 2 
We have the following result whose proof can be found in Appendix D. 

Lemma 6.7. 

\Roo(t)\ < Roo(0), \R'Ut)\ < |iC(0)|, Vt > 0, 



Rooit), i4(t), R'U{t) = O (jj as t 



oo. 



Roo (t) 



(A^-OA 2 -A (^ 00 ) i >0, Vi>0. 
^ 2 + ^ 4 + 0(t 6 ), J R' 0O (t) = -^ + ^ 3 + 0(t 5 ), ast 



(6.19a) 
(6.19b) 
(6.19c) 

(6.19d) 

□ 



We set 

S(t) := max(i, 1), t > 0. 
We find it convenient to introduce new functions 



R,{o) 

Using these notations we can rewrite (6. 19c) as 



j^RAt), Hv {t) = ^ ) K{t) = - J ^ 7) 



Kit)- 



(i-gL)-^(g / oc ) 2 >o, Vi>0. 



A 2 



:r,(t) 



The equalities (6. 19d) imply that 



Then 



where 



and 



fu(t) 



v(t) 



A 2 H 



4375 



t A + o(t 6 ), v|t| < i. 



x Q v (t) x - p 2 y + p u arcsinp^ 



(1 - G u {t?) - ^(G' u (t)) 2 



(1 - G v (t? fl 2 

K(X^?-Rl) + (K?R, . -Hu(l-Gl) + ^(G' u rG u 



(6.20) 



(6.21) 



yv (Ag^-^^-AoC^) 2 
Lemma 6.8. Let k G (0, 1). Then 

e u (t) = 0(t), V0 < t < v~ K , 
where the constant implied by O -symbol is independent ofv and t G [0, v~ K ), but it could depend on 
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Proof. Observe that for t G [0, v~ K ] we have 

G ^= i -^ 2+ °( t4 )- a ^=-W> t+o{t3) 

so that 

-3 



and 



(l-G,(t) 2 )-^-(G'J = 0(t 4 ) 
so that Q u (t) = O(i). □ 
Lemma 6.9. Le? «e (0,1), 77ze?2 



e,(t) - CooCt) + O ( J 7( ^3/2 + uS ft\l)3/2 + ^( t)7 (t)3/2 + ^3/2 ) ' (6 " 23) 
»7(t) := (1 - G 2 ) - ^(G'J 2 W 7 (t) = 1 - G^tf. 



Proof. Observe that 



so that 

l-G OQ (t) 2 (t) = (1-Goo(0 2 ) (l + O 



i// y 00 \vS(t) 

i 



^(*)7(<) 



For t > v K we have 



i/<S(*)7(t) 
Hence 



/ l 

O I — j— - I = o(l) uniformly in f > v K as — > oo. 



[I - G 2 )^ 3/2 = (1 - G^)^ 2 (l + O (^^y) ) , 70 := 1 - Gl(t). (6.25) 



Next observe that 

3z> ~ 3 ' ' " _ 5i^ 5 5 



Ao(-) = ^i^ = U-- 1 + 0(,- 2 ), A 2 H = ^±11 = 1 + ^ + 0(0 (6.26) 



and 

1 i ,,-1 i /0/,,-2 



AoH _ l + ^-i + Q(^) _ 5 _ 10 , 2 
AzH ± + ^ + 0(^-2) 3 3 " )■ 
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Using (6.24) and the above estimate we deduce 

A o^) (r> \ 2 = n _ r 2 \ - 



(1 " GSf) - Wk(G' v f = (1 - Gl) - ^(O 2 +^(0 2 + O ( -L- )+ O ( 



(1 _ _^).( 1 + (^ + ^ + _l_)) 



(6.27) 

Using (6.25) we deduce (6.22). The estimate (6.23) follows (6.22) by invoking the definitions of j(t) 
and 77(f). □ 

Lemma 6.10. Let k £ (0, j). Then 

( 1 (G' ) 2 1 1 1 \ 

Pv = Pco + \Mt) + JZffis® + + rffflv® + ^WW) ) ' yt>u ~ K > (6 " 28) 

where the constant implied by O -symbol is independent of v and t > u~ K , but it could depend on k. 
Proof. The estimates (6.19b), (6.19d) and (6.21) imply that for t > v~ K we have 

M°2! + I + 1 = or 1 ) = 

uniformly in t > u~ K . We conclude from (6.27) that 



( (1 - cm 2 ) - j^j^y')' 1 = ( f 1 - °y - ^(<? 



0^/ x2 

00 / 



1/77 (t) vS(t)r){t) v 2 rj(t) 



(6.29) 



Since ifj, = i?oo +0{v ) we deduce 



-H„(l - G v f + M4(G'J 2 G, = -^(1 - G^) 2 + t^(G" 00 ) 2 G 00 + O f- 

A2(fJ A2 \^ 



Recalling that 

u n r>2^ j 

A 2 ( 



Pp 



H v {l-Gl) + ^{G' V ) 2 G V 



^(rJ \ 2 n - n(x- l \ 

Consider the function 



and -#oo(l - Goo) 2 + T a (G' 00 ) 2 Goo = 0(6^) we see that (6.28) follows from (eq: est22). □ 



A(u) = yl — u 2 + u arcsin u, |u| < 1. 

Observe that 



(L4 

— = 0(1), Vlttl < 1, (6.30a) 
du 

dA 

— = arcsin u = 0(u) as u — > 0. (6.30b) 

du 



Now fix an exponent k G (0, 4). We discuss separately two cases. 
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1. If v k < t < ttu. Then in this range we have 

Poc = 0(t~ 1 ), Pv = Po0 + 0(l), ^y = 0(l), 

and using (6.30b) and (6.28) we deduce 

1 

2 



A{ Pu ) = A{ Poo ) + A'( Poo )( Pu - poo) + 0( (p u - poo) 2 ) = A( Poo ) + O + 
2. v~ K < t < v K . The equality (6.21) shows that in this range we have 



. (6.31) 



v (t) " " 8{t) 

so that and (6.28) implies that 

Pu~P^ = {-^ 

Using (6.30a) we deduce 

A( Pu ) - A( Poo ) = O (-^ + -— ^ + — L_ ) . (6.32) 

Set 

A v := e u A( Pu ) - eooA(poo), q v :=- 



Then using (6.26) we deduce that 

6 



5 



v- 1 + 0{v~ 2 ), q v = goo (1 + 2u~ x + 0(z>- 2 ) ) . 



Then 



fu(t) - q v J - ( /«>(*) - Qoo J = 

= ^(eoo^oo(Poo) - 1 + A„) - 5oo(eoo^4(Poo) - l) 
= (g V - <?oo)( eoo-4(Poo) - 1) + qoo^v 

To prove (6.1) we need to prove the following equality. 

\Qv 9oo) / ( Sqo 

(t)A( Poo (t)) - l)dt, qoo / A v (t)dt = o(l) as i/ ->■ oo. (6.33) 

We can dispense easily of the first integral above since C 00 (t) J 4(p 00 (i)) — 1 is absolutely integrable 
on [0, oo) and q v — q^ = 0(v~ l ). 

The second integral requires a bit of work. More precisely, we will show the following result. 

Lemma 6.11. I/O < k < g, then 

/ A u (t)dt, / A u (t)dt, / A v (t)dt = o(l) <w f ->■ oo. (6.34) 

Proof. We will discuss each of the three cases separately. 

1. < t < v~ K . The easiest way to prove that A v {t)dt — > is to show that 

e v (t)A( Pv (t)) = o(i), o<* 

This follows using Lemma 6.8 and observing that the function A(u) is bounded. 
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2. v K < t < u K . In this range we have 

l - = 0(v^% _L = (^), J^ = 0(1), G' 00 (t) = 0(l). 
Using (6.23), (6.32) we deduce 

Qu(t) = eoo(i) + O (-^) , MPu) = A( Poo ) + O [y^j ■ 



Hence 



A u = o(^^j and J^A u (t)dt = o(j-^=o(l) 



3. v K < t < wis. For these values of t we have 



Using (6.23) and (6.31) we deduce 

e„(t) = eoo(t) + o f 1 + j , A( Pu ) = A( Poo ) + o + 1 

so that 

M*) = O f i + ^ V /™ A„(t) <ft = O f - + = o(l) 



v v 



□ 



The fact that defined as in (6.1) is positive follows by arguing exactly as in [19, §3.2]. This 
completes the proof of Theorem 6.1. □ 

Remark 6.12. The proof of Lemma 6.11 shows that for any e > we have 

var{Z v ) = udoo + 0(v e ) as v — > oo. (6.35) 
Numerical experiments suggest that 5^ « 0.35. □ 

Appendix A. Some elementary integrals 
Suppose W is an oriented Euclidean vector space equipped with an orthogonal decomposition 
W = Wo®Wi, dimWi = ni, i = 0, 1, dim W = n = n + n\ . 

For any w 6 W we denote by wi its orthogonal projection on Wi, i = 0,1, so that w = wq + w\. 

We set ri{w) := \wi\, i = 0,1. 

Lemma A.l. Let ifi : Wi — >■ M, i = 0, 1, be locally integrable functions, such that (p$ is positively 
homogeneous of degree ko > 0, and set (p(wo,W\) := (po(wo)ifi(wi). Then 



<p(w)\dS(w)\ = <p (wo)\dS(w )\x ^(w^l-rO^^ldViw^. 

<S(W) JS(W ) Jb{W!) 

(A.l) 

In particular, if ip$ = 1, then 

[ ip(w) \dS(w)\ = o- no _! f ^(w^il-rf^ldViw^l (A.2) 
JS(W) JB(Wi) 
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Proof. The key trick behind the equality (A.2) is the co-area formula. Denote by ir the orthogonal 
projection onto W\, i.e., 

7T : W — > W\, Wo + Wl !->■ W\. 

This induces a smooth map it : S(W) — > B(W\). We denote by Y, Wl the fiber of this map over 
wi. We observe that T, Wl is the sphere in Wo of radius (1 — r 2 ) 1 / 2 and center 0. Denote by 
J w : S(W) — > M the relative Jacobian of the map it defined as in [22, §5.1.1]. 

Fix a point w = (wq, W\) 6 S,^. Next, choose an orthonormal basis ei, . . . , e rei of Wi such 
that ei = —wi . 

i rj 1 

The orthogonal complement N w T, Wl of T w Y< Wl in ^^(W) consists of vectors that are orthogo- 
nal on Wo and on the unit vector w = wo +W\. We deduce that the collection 

fi = -rfwo + rlwx, f 2 = e 2 ,...,f ni = e ni , 
is an orthogonal basis of N w Y* Wl . Note that 

I n I 4 2 , 4 2 2 2 
l/ll = r l r + r r l = r l r 2- 

We obtain an orthonormal basis by replacing / 1 with the vector 

l/il n 

The orthogonal projection onto W\ of the orthonormal basis / 1; / 2 , . . . , / ni is the orthogonal basis 

r ei,e 2 , . . . ,e ni , 

whose determinant is ro- This shows that 

J w ( w ) = ro {w) = (1 - niw) 2 ) 1 / 2 , Vw e S{W). 
The coarea formula [22, Thm. 5.3.9] implies that 



<p(w)\dS(w)\= f if - -(p(w ,wi) \dSiw! 



S(W) Jb{Wi) \Jx wl Jtt{wo,Wi 

<Po(wo)\dS(w )\) [ ^{w^r^-^dViw^ 
S(W ) J JB{Wx) 

^o(^o)l^o)l) / Mwi)(l-r 2 1 ) tQ± ^ 1 \dV(w 1 )\. 
S(W ) J JB{W\) 



□ 



Suppose that L is a Euclidean vector space of dimension £, and Q : L — > R is a continuous, 
positively homogeneous function of degree k > 0. For any positive integer n we set 

In(Q):= [ \Q(x)\(l-\x\ 2 ) n / 2 \dV(x)\, J n (Q):= [ \Q(x)\\dS(x)\, 
Jb(l) Js(l) 

where S(L) denotes the unit sphere in L centered at the origin, and B(L) denotes the unit ball in L 
centered at the origin. 



Lemma A.2. 

Y(i±h)Y(z + i) 
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Proof. We have 



In(Q) = UQ) C r i+k -\l-r 2 ) n l 2 dr= l -J n {Q) f s^-\l - s) n / 2 ds 
Jo 1 Jo 

where B(p, q) denotes the Eulerian integral 

T(p)Y{q) 



B(p,q) = f (l-sf-V^ds 
Jo 



T(p + q) ' 

□ 



The integrals of homogeneous functions over the unit sphere can be reduced to gaussian integrals 
of these polynomials. More precisely we have the following result. 

Lemma A. 3. Suppose that W is an Euclidean space of dimension N and f : W —> R is a locally 
integrable positively homogeneous function of degree £ > 0. Then 

I f(x) \dS(x)\ = —^j- I e-\^f{x) \dV(x)\, (A.4) 
Js(w) 1 {— Jw 

— — / f(x) \dS(x)\ = -^L / f( x )L^- \dV(x)\. (A.5) 

area S(W) J s{w) r(^±£) Jw vrT 



or equivalently, 



Proof. We have 

/ e-W 2 f(x) \dV(x)\ = Hi [ f(x) \dS(x)\) e~ T \ N ^dr 
Jw Jo \JS(W) / 



/ f{x) \dS{x)\ \ r e-'s^- 1 ^ = [ fix) \dS(x) 

JS{W) J Jo 11 JS(W) 



□ 



Proposition A.4. Suppose a and b are nonnegative real numbers such that a > b. Then 

I(a,b):= [ e-^^ {ax2+ay2 - 2bxy) -^ 2 \xy - z 2 \\dxdydz\ 
Jr 3 

, / r 27r 2c 3 / 2 \ 

= V2vr(«2-62) U 27ra + 27rj , 

where c{9) := (a — b cos 28). 

Proof. Let {i,j,k} be the canonical orthonormal basis of M. 3 . Define a new orthonormal basis 
e\, e 2 , e3 of ]R 3 by setting 

ei = — ^(i + i); e 2 = -j=(i-j), e 3 = k. 

If we let (u, v, w) denote the coordinates with respect to this new orthonormal frame, then from the 
equality 

ue± + ve 2 + we3 = xi + yj + zk 
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we deduce 

z = w, x = -j=(u + v), y = -j=(u-v) 



xy = -^w 2 ~ ^ 2 )' ax2 ~ ^bxy = (a — b)u 2 + (a + b)v 2 . 



We deduce that 



I = - e 2 ( a2 - i ' 2 ) 2 | it - v - 2wT| |cfat(fodi(;| 

We now make the change in variables 

u = ■\/2(a + b)u, v = \/2{a — b)v, w = V2w, 



to deduce 



/ = 2^2(a 2 - b 2 ) [ e -(« 2 +^ 2 +«' 2 ) | ( a + b)u 2 + (a - b)v 2 - 2w 2 \ \dudvdw\ . 



We now change to cylindrical coordinates, 

w = w, u = rcos9, v = rsin9, 

so that 

v 2 + v 2 + w 2 = r 2 + w 2 , 
(a + b)u 2 + (a - b)v 2 - 2w 2 = r 2 ((a + b) cos 2 9 + (a - b) sin 2 6 ) - 2w 2 
= r 2 (a - bcos26) - 2w 2 . 

We have 

i-2-k poo roo 

II 



f'ZTT f'OG f'OO 

/ del dw j e~ r2 1 r 2 (a -b cos 26>) - 2u> 2 | rdr 

JO J-OO JO 

- d6 e~ w ' 2 dw e- s \s(a-bcos29)-2w 2 \ds 
2 Jo J-oo Jo 



At this point we observe that for any c, d > we have 

roo d 
/ e~ x \cx — d\dx = 2ce~~ + d — c. 
Jo 

Hence, if we set c = c(9) = (a — b cos 29) we deduce 



1 r2n roc 2 

Ii = - d9 e- w2 (2ce-^ +2w 2 - c)dw . 

2 Jn J-00 



=:J(c) 

We have 



J(c) = 2c I e- £ ? w2 dw + 2 [ e- w2 w 2 dw-c [ e~ w2 dw 
Jr Jr Jr 

7T 1 / 2 + 21X3/2) - C7T 1 / 2 = 7T 1 / 2 ( - 2c3/ ' /9 - C + 1 



(c + 2)V2 w , ^( c+ 2)1/2 

We deduce that 

2c 3 / 2 



/*27T /*27T / 

/ = ^(a 2 - b) 2 J J(c)d9 = v / 2TT(a 2 - b) 2 J ( 



(c + 2)V2 



c + 1 A 



The conclusion of the proposition follows by observing that c(9)d9 = a. □ 
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Appendix B. Basic facts about spherical harmonics 

We survey here a few classical facts about spherical harmonics that we needed in the main body of 
the paper. For proofs and more details we refer to our main source, [25]. 

We denote by "K n d the space of homogeneous, harmonic polynomials of degree n in d variables. 
We regard such polynomials as functions on M. d , and we denote by ^ n ^d the subspace of C°°(S d ~ 1 ) 
spanned by the restrictions of these polynomials to the unit sphere. We have 

,. nr . . fd + n—\\ fd + n — 3 N 

dim5i nid = dimy n>rf = M(n,d) 



Observe that 



n J \ n — 2 

2n + d-2(n + d-2\ n d - 2 

= : ~ 2— — as n — > oo. 

n + d-2 \ d-2 J (d-2)\ 

M(0,d) = l, M(l,d) = d, M(2,d) = ( d+ 2 l \ - I. (B.l) 



The space y n ^ is the eigenspace of the Laplace operator on S d 1 corresponding to the eigenvalue 

X n (d) =n(n + d- 2). 

We want to describe an inductive construction of an orthonormal basis of Vin,d- We start with the 
case d = 2. For any m G Z, we set 

, n . fcos(m0), m<0 f(27r) 1/2 , m = _ 1 

<Pm(V) = < . , n , t m = \\ip m \\ L 2 = < I/O n ) *m = — ¥>m- 

Then So,2 = {^o} is an orthonormal basis of ^0,2, while 23 n; 2 = {3>-n> ^n} is an orthonormal basis 
of V n ,2, n > 0. 

Assuming now that we have produced orthonormal bases < B n ^_ 1 of all the spaces )$ n d-x, we 
indicate how to produce orthonormal bases in the harmonic spaces y n ^- This requires the introduction 
of the Legendre polynomials and their associated functions. 

The Legendre polynomial P n ^d{i) of degree n and order d is given by the Rodriguez formula 

where R n (d) is the Rodriguez constant 



rc^=i) 1 

R n (d) = 2~" - 1 2 / 1x = 2~" 



r( n+^i) [n+^]„' 

where we recall that [x]k := x(x — 1) ■ ■ ■ (x — k + 1). Equivalently, they can be defined recursively 
via the relations 

P , d (t) = l, Pi,d( t )= t i 
(n + d- 2)P n+M (t) - (2n + d - 2)tP n>d (t) + nP n _ M (i) = 0, n > 0. 
In particular, this shows that 

P 2>d (t) = ^(dt 2 -l). 
The Legendre polynomials are normalized by the equality 

P n Jl) = l, Vd>2, n>0. 
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More generally, for any n > 0, d > 2, and any < j < n, we have 



^2(1) = i-V n Rn{d) 



n + j\\ iff (1 ~ t) n 
J J\ (l-t) 



Di(l+t) 



d-3 
2 



d-3 
2 



D f := — 



which implies 



2 n - j R n (d) 



n + j 
j 



d-3 



n 



(1+*) 

d-3 



t=i 



dt' 



n + 



n + j 
j 



n + 



d-3 



(B.3) 



For any d > 3, n > and < j < n, we define the normalized associated Legendre functions 
where 



a 



[n + d- 3] d -3 



(2n + d- 2) 



r( d_i ) ^ 2 <*-> + d + i-3] 2i+d _ 



1/2 



When d = 3, the above formulae take the form 

Pn,j,3{t) = 



{n + \){n-j)\ 



(B.4) 



(B.5) 



V (n + j)! 

For any < j < n, and any d > 2 we define a linear map 

^n,j,d '• ^j,d—l ~ ^ ^n,di Y I— > 7 n j id [Y], 

7 nJ4 [Y](x) =Pi/x d ) .yL^-x'J, VxeS d -\x' = (x 1 ,...,x d „ 1 )^0. 

Note that for x = (x', x d ) G S"^ 1 we have 

\\x'\\ = (1 - x,) 1 / 2 and Pl d {x d ) = C nd , d (l - x d y' 2 P^ d (x d ) = C n>j , d \\x'\\ j P^(x d ), 
so that 

7 nJ4 [Y](x) = C nJ4 P^ d (x d )Y(x'), Vx = (x',x d ) e S d ~\ 

where Y denotes the extension of Y as a homogeneous polynomial of degree j in (d — l)-variables. 
The sets T nj - dpBj,<f-i], < j < n are disjoint, and their union is an orthonormal basis of ^ n>d that 
we denote by B n d . 

_ i 

The space y ,d consists only of constant functions and "Bod = { a d ~\ }• The orthonormal basis 
of obtained via the above inductive process is 



^>\,d = { C Xi, 1 < i < d] = {cr d 3 2 C lt0>d Xi; l<i<d}. 
The orthonormal basis B2 d of ^2 d is 



Ci(dx 2 



1 < i < <i, C2XiXj, I < i < j < d, 



(B.6) 
(B.7) 



where r 2 = x 2 + • • • + a^, and the positive constants Co, Ci, C2 are found from the equalities 

Cl [ x\ \dS(x)\ = C\ [ (d 2 x\ - 2dx\ + 1) \dS(x)\ = C\ [ x\xl\dS(x)\ = 1, 
Js d - 1 Js d - 1 Js 6 -- 1 
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aided by the classical identities, [27, Lemma 9.3.10], 



x 2 ^---xf d \dS{x)\ = V 2 J 2h+d , 2 , h = h 1 + --- + h d . (B.8) 



2/ii 2h d I in/ \i _ ^ V 2 

Ju 

S d-1 



Appendix C. Invariant integrals over the space of symmetric matrices 
In the main body of the paper we encountered many integrals of the form 

| det v4| \drf(A)\, 



'Sym N 

where Sym^ is the space of symmetric N x N matrices, and 7 is a Gaussian probability measure on 
Symjy. In this appendix, we want show that in certain cases we can reduce this integral to an integral 
over a space of much lower dimension using a basic trick in random matrix theory. We set 

D N :=dimSym N = (^^y Sym% := { S € Sym N ; tr,S = 0}. 

Note first that we have a canonical 0(iV)-invariant metric on Sym^ with norm | — |* given by 

\A\, :={trA 2 ) 1/2 . 

Using the canonical basis of M. k we can describe each A £ Sym^y as a linear combination 

A = a ijHij, 

where Hij is the symmetric k x k matrix whose and entries are 1, while the remaining 
entries are 0. With respect to the coordinates (a^- ) we have 

g* = J2 da 2 + da%. 

i i<j 

The collection (-ffy)i<i<j<Ar is an orthonormal basis with respect to the metric g*. The volume 
density determined by the metric g* has the description 



\dV* 



~[dai 

i<3 



Via the metric on M. N we can identify Sym^ with the vector space of homogeneous polynomials of 
degree 2 in iV- variables. More precisely, to such a polynomial P we associate the matrix Hess(P), 
the Hessian of P at the origin. The subspace Sym^ corresponds to the space !K2,7V of homogeneous, 
harmonic polynomials of degree 2 on R . 

The orthogonal group 0(N) acts by conjugation on Sym^y, and Sym^y decomposes into irre- 
ducible components 

Sym^v = R(1jv) © Sym^, 
where R(ljv) denotes the one-dimensional space spanned by the identity matrix Ijy. 

We fix an 0(iV)-mvariant metric on Symjy. The irreducibility of Sym^y implies that such a 
metric is uniquely determined by two constants a, b > so that the collection 

al N , 6Hess(y), Y G S 2 ,tv 

is an orthonormal basis. We denote by | — | a ,& the norm of this metric. We want to express |-A| a ,& in 
terms of tr A 2 and (tr A) 2 . 
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Note first that 



\ 1 »\l fi = i i = 1 ^i\iN\l. 



The irreducibility of Sym° implies that there exists a universal constant R = > such that for 
any homogeneous harmonic polynomial P of degree 2 in N variables we have 



If we take P = x\x%, we deduce 



Hess(P)| 2 = R 2 [ P(x) 2 \dS( 
x\xl\dS{x% 



2 = R? 



S N-1 



and using (B.8), we deduce 



We see that for any P £ IK 2 



R 2 



N 



r( /v±4 } 4 r(^) 



r(3/2) 2 r(i/2)^- 2 
l 



N 
7T 2 



(C.l) 



\ Y \a,b ~ ^1^11x2(5^-1) - b2R2 



Hess(r)| 2 . 



In particular, we deduce that 



1 h -I 

N' b *~ R- 



In general, if A G Sym^, then we have a decomposition 

A = — (ixA)l N + (A- -i(t 



that is orthogonal with respect to both 



and 

2 



^(trA)ljv 



a.b 



N 

\ a fi. We deduce 
+ (A-l(trA)ljv; 



( tr ^) 2 + 7^ tr ( ^ - |(tr A)A + ^(tr A) 2 1^) 



iV 2 a 2 



1 



6 2 i? 2 
1 



P^) (trj4)2 + P^ tr>l2 - 



(C.2) 



::/3 



Note that the quantities a, (3 depend on a, b and the dimension N. 

If | dV* | denotes the volume density determined by the metric | — | * and | dV a ^ | denotes the volume 
density associated to the metric | — \ a ^, then we have 

\dV a>b \ = C N (a,b)\dV*\, C N (a,b) := • (C3) 

Suppose now that / : Sym^y — > E is a continuous 0( 7V)-invariant function that is homogeneous of 
degree £ > 0. We want to find a simpler expression for the integral 



Sym^ 



■^■»f(A)\dV atb (A)\. 



This can be reduced to a situation frequently encountered in random matrix theory. We have 



J a> b{f) = C N {a,b) 



Sym f 



-^f(A)\dV*(A)\ = C N (a,b) 



Sym f 



- atrA2 -^ A ^f(A)\dV*(A)\ 
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Via 



iSym w VV2a / fO-TT^ 2 



=:*a,js(S) 

Observe that the function f? i-> <& a ^(B) is also O(N) invariant. We denote by Djy C Sym^ the 
subspace consisting of diagonal matrices. We identify D n with in the obvious fashion. Using [2, 
Prop. 4.1.1] or [12, Thm. 2.50] we deduce that 

75^ 1^*1 = ^ / ®a,p(B) |A(S)| ^TOI, 

's ymiV (2tt)— Zjv ^iv (2vr) — 

where 

• A (a;) is the discriminant A(xi, . . . , Xjv) = IIk^x/vC 3 ^ ~~ x j) 

• The constant Zjv is given by the integral 



/ 

J Si 



- [ \A(B)\ ■ e -^L\dV{B)\ = (2tt)-t / lA^lf-^ \dV(x)\. 



(2tt) 2 ^ (2tt) : 



* ^ r(i + 1) 



n 

j=0 ^2 



Putting together all of the above, we deduce 

D N 

JaAf) = ^ CN ^ b) [ e-^-^-^ff (*) \A( X )\ \dV( X )\. (C.4) 
Z*m-f J^ N \V2aJ 



In particular, we have 



e -\ A ^\det{A)\\dV a ^A)\ 

Symjv 



7T 2 C/yfa, b) f M 2 §_(s^ N -r .12 



N 

Zn(2q)~ a~ 



N 



l\xi\-\A(x)\\dV(x)\, (C.5) 



i=l 

where a, ft are defined by (C.2) and Cat (a, 5) by (C.3). 

Let us point out that, up to a universal multiplicative constant, the measure e~2 trA \dV\*(A) is 
the probability distribution of the real gaussian ensemble, [2, 12]. As explained in [12, Chap.3], 
the multidimensional integral (C.5) can be reduced to computations of 1-dimensional integrals in the 
special case when ft = 0, i.e., ka 2 = b 2 R 2 . As explained in [17, §1.5], [18], the case ft < can 
be reduced to computations of 1-point correlations of the Gaussian ensemble of (k + 1) x (k + 1)- 
matrices. In turn, these can be reduced to computations of 1-dimensional integrals [12, §4.4], [17, 
Chap. 6], [23, Chap. 7]. 

Appendix D. Some elementary estimates 
Proof of Lemma 6.2. Consider the complex valued random process 

X — a tmt 

3^v\t) • — i / TfiZm,& u , z m — c m id m . 

m =l 
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The covariance function of this process is 

1 v 

= E(%(t)%{0) ) = — V 

m=l 

Observe that Re 3 U = <Ev Note that the spectral measure of the process "J v is 

1 " 

d&v = — o 

r7/o / 



m 2 5i 



m=l 



where St denotes the Dirac measure on R concentrated at to- We form the covariance matrix of the 
gaussian vector valued random variable 



9v(0) 
3v(*) 



►X„(t) 



B„(0) 

D^(t) 3^,(0) 



31,(0) 

-Kit) 



Kit) 

-<(0) 

-<(t) 
-K(o) 



Observe that Re X(i) = Re E(i). If we let 



Then, as in [9, Eq. (10.6.1)] we have 

1 



(X u z, z) 



3 Z_^i 

m=l 



m 



"0 

III 



u + v e~ 



+ 



im 



We see that 

{X u z, z) = 0- 



u + foe * 



H I ui + v±e <> 

v 



u\ + v\e~ 



0, Vm = 1, 



(D.l) 



We see that if the linear system (D. 1) has a nontrivial solution z then the complex 4x4 matrix 



i c 

i c 2 

i c 3 

l c 4 



Mt) ■ 

must be singular, i.e., det A u (t) = 0. We have 
det A v {t) = det 



1 c 

2 2C 2 

3 3C 3 

4 4C 4 



e " 



l c i c 

o c 2 -C 1 2C 2 -z 
C 3 -C 2 3C 3 -z 
— z 



C 2 det 



c 4 

c- 
c 2 - 
c 3 - 



2 



3 AC 4 



1 1 2C - 

-1 2 3C 2 
■1 3 4C 3 



1 
1 
1 



C 2 det 



C 2 det 



C 3 det 



C- 1 i i 

C 2 - 1 2 2C 
C 3 -l 3 3C 2 



Cdet 



c- 

C 2 

c 3 



" 1 


1 1 


1 







C-i i 


2C- 


1 





C 2 - 1 2 


3C 2 " 


1 





C 3 - 1 3 


4C 3 - 


1 


c- 


-liC" 






c 2 - 


- 1 2 2C 2 






c 3 - 


- 1 3 3C 3 








1 








2 2C 2 - 2 








3 3C 3 - 3 
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C 3 det 



C - 1 1 
C 2 - 2C + 1 2C 2 - 2 
C 3 -3C + 2 3C 3 -3 



C 3 (C"l)det 



C det 



1 1 

= C 3 (C"l) 3 det 1 2(C + 1) 

_C + 1) 3(C 2 + C + 1) 

= -C 3 (C-i) 3 (C 2 -C-i)- 

Since |C| = 1, we that det A v (t) = if and only if t G 2-kvI,. 



C- 1 1 

(C-i) 2 o 2 (c-i)(C + i) 

C 3 - 3C + 2 3C 3 - 3 

C - 1 1 

(C-l) 2(C + 1) 

(C-i) 2 (C + i) o 3 (C-i)(C 2 + C + i) 



C 3 (C-i) 3 ( 3(C 2 + C + i) - 2(C 2 + 2C + 1) ) 



□ 



Proof of Lemma 6.6. Recall that 9 :- 



2v 



./(*) 



t r+l 



2 sin 



D„, r (t) = r! 



(v+l)t 
2u 



r+l 



i(f + r)t 

e 2" 



. By (6.10) we have 

r +1 ) 



3=1 



V3 



r+l-j 



f(0)J 



Using (6.13) we deduce that 

t r+l 



D v , r {t)) ~ -I r {t) 



< 2rH r+1 



sm 



f(oy +1 ' 



vi \f(9) 



(D.2) 



+ r! ^ t j 

3 = 1 

In the sequel we will use Landau's symbol O. These implied constants will be independent of v. 
Also we will denote by the same symbol C r constants independent of u put possibly dependent on r. 
Throughout we assume < t < ttu. Then < 9 < ^ and for < j < r we have 



jo 



1 + 0(9), 



(T) 



[ - 

V 



sm 



9(cos9 + isin#) 



Hence 



f(0) 



sm 



2u 



sm I 



e ire - sin 



while for any 1 < j < r we have 

(v+l 



v j \fiP). 
Using (D.3) and (D.4) in (D.2) we deduce that 



mi 



%e \ r+i-j 



> + !)* 
2v 

1 + 0(9). 



o 



1 



+ 



(D.3) 



(D.4) 



t r+1 1 

—p^D^t)) - -:I r (t) 



% V 



3=1 



3=1 
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Hence 



1 



Ir(t) 



o 



1 (1 - t r 

~vF (i-t) 



□ 



Proof of Lemma 6.7. We have 



R&Ht) = ± 



1 



T k+2 u(r)dT, u{t) 



Note that 



i 



j sinr, k £ 1 + 22 
1 cos t, A; £ 2Z. 

r fc+2 dr. 



" jfc + 3 t k + 3 

The inequality (6.19a) now follows from the inequality \u(t)\ < 1, Vr. 

For any positive integer r we denote by j r the r-th jet at of a one- variable function. We can 
rewrite (6.13) as follows: 



1 



I r {t) = rl 2 sin 



t\ a 
e 2 — xe 



r! ( 2 sin 



i \ it 



e 2 — ^e 



•jV^-i) 



r! I 2 sin I - I e~ + ie** - ie l * • j r ( e~ l * ) 



Hence 



rl I — z( e 2 — e 2 



1 



)e 2 + ie rf - ie 1 * • j r (e~ 



-it 



Re -J r (t) =Im e Jt -j r e 



and 



ir! 1 - e • j r ( e 



I r (t) = 0(t _1 ) t ->■ 00. 



This proves (6.19b). 
The spectral measure 



1 y 



m 2 bm. 



m=l 



of the process "J v converges weakly as v — > 00 to the measure 

— X[<d,i}t 2 dt, 



da r . 



7T 



where X[o,i] denotes the characteristic function of [0, 1]. Indeed, an argument identical to the one used 
in the proof of Lemma 6.4 shows that for every continous bounded function / : R — > R we have 



lim 

v— >oo 



f(t)d<j v (t) = / /(t)dff, 



The complex valued stationary Gaussian process 1^ on E with spectral measure (fcoo has covariance 
function 



- / t 2 e u dt. 

7T 



Note that Re^oo = R^. The results in [9, §10.6] show that the covariance matrix 

^oo(O) 3?oc(t) 3C(0) ^o(t) 
Koo(t) ^cc(O) -iK'ooCt) -^oo(O) 

^oo(O) *<o(t) -3^(0) 



lim Xjy, 

1/— >oo 
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is nondegenerate. The equality det Re Xoo(t) ^ 0, Vt £ R implies as in Remark 6.3 that Hoc(t) ^ 0, 
| Poo (01 < 1, Vi € R, where 

_ (Ap - -^00)^2 - A (^) 2 -fi"4(Ap ~ j^go) + ( r 'oo) 2r oq 

Ag-l^ ' (Ag-^)A 2 -Ao(^) 2 • 

This proves (6. 19c). The equality (6. 19d) follows from the Taylor expansion of R^. □ 
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